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Course summary 


► Reservoir simulation workshop Introduction 

► Reservoir simulation workshop Simulator presentation 

► Dynamic reservoir simulation Data review: gridding 

► Dynamic reservoir simulation Data review: PVT 

► Dynamic reservoir simulation Data review: Initialization 

► Dynamic reservoir simulation Data review: Aquifers 
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► Dynamic reservoir simulation Data review: Flow description 

► Dynamic reservoir simulation Data review: Well representation 

► Reservoir simulation workshop History match 

► Reservoir simulation workshop Production forecast 
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An IFP Training Course for REPSOL 

Reservoir simulation workshop 
Introduction 

Instructor: Maria AGUILERA 

/FP7 'raining 


Course outline 


► Introduction 

► Simulator presentation 

► Data review 

• Reservoir description - Gridding 

• Fluid description 

• Initialization 

• Aquifers 

• Flow description 

• Well representation 

► History matching 

► Production forecasts 
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Objectives 




Objectives 



► Introduce you to the basic concepts of modelling 

► Introduce you to the numerical reservoir modelling capabilities, and major 
features 

► Illustrate with some results, with examples and exercises 

► By the end of this lecture, you should be aware that: 

• Simulation is a useful reservoir engineering tool 

• A model must not be used as a black box 

• The output results depend upon the quality of the input data 

• What goes in what comes out!!! 
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► What is Dynamic Reservoir Simulation? 

• Reservoir simulation is an area of reservoir engineering in which computer models are 
used to predict the flow of fluids (typically, oil, water, and gas) through porous media, 
before deciding on a development plan and investment strategy 

• Simulation is a reservoir engineering tool; it is the most powerful tool in the reservoir 
modeling tool box. 

• Reservoir simulation associated with material balance study are the basic techniques an 
engineer must be familiar with. 
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Dynamic reservoir simulation 


► Early in a fields life, it is essential to be able to evaluate how much 
Hydrocarbons will be produced through time: 

• How much hydrocarbon is in the reservoir? 

• How much of that hydrocarbon can potentially be recovered? 

• How quickly can the recoverable hydrocarbons be produced? 

• How will the reservoir perform under various development scenarios? 

► A Dynamic Reservoir Model is a useful tool that may help to answer those 
questions. 
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Why reservoir simulation... 


► Main drivers: 

• To overcome simplifying assumptions (used in hand calculations) 

• To consider realistic problems (development plan, production history) 

• To take into account all the available data and to better understand all the interactions 
between them for a given scenario 

• To perform sensitivity to unknown parameters 

► Main use: 

• Help to decision, in particular at the end of appraisal phase and during pre-development 
screening 

• Optimize the production profile and the final recovery after production start-up by 
improving the reservoir description 

• Establish of reliable production forecasts 
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Why reservoir simulation... 



► Appraisal 

• To identify dynamic incertitude 

• To define appraisal needs (contacts, faults, facies variations) 

► Non producing reservoir: Field development plans 

• To identify efficient recovery mechanisms 

• To define the well number, their location and their architecture 

• To give production profiles to decision makers 

• To appreciate the risks associated to the project 

• To optimize Capex 

► Producing reservoir or Mature fields 

• To improve the reservoir characterization through history matching 

• To give production profiles (oil, gas and water flow rates, pressure, fluid composition, ...) associated to 
existing wells 

• To define completion changes or infill drilling 

• To optimize oil production and recovery 

• Optimize field capacities (treatment capacities, compressors, ...) 

► In all cases 

• To establish reliable production forecasts 
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Reservoir simulation... 
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Reservoir simulation model 
Data set workflow 



► Reservoir characterization 

• Depositional environment, structure, rock and fluid properties 


► Geological model 

• Fine grid model filled up by automatic distribution methods; calculation of hydrocarbons 
in place volumes 


► Upscaling process 

• Average small cells properties to fill up larger cells 


► Reservoir simulation model 

• Coarse grid model for history match and prediction runs 
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Reservoir modeling workflow 
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Upscaling 


► Geological model contains a very 
large number of grid blocks. This 
level of resolution is generally 
incompatible with the computing 
capabilities of numerical flow 
simulators. 


► To do so, fine grid blocks are 
grouped in aggregates: coarse grid 
blocks. The basic problem is to 
determine the equivalent properties 
of the coarse gridblocks. 



Fine grid 
(geological model) 





Coarse grid 

(reservoir model) 
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Reservoir simulation model 


Fluid Characterization 
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Reservoir simulation models 


Reservoir simulation model 


Numerical simulator 

+ 

Set of data 

► Modeling: building the model (define parameters through integrated 
reservoir characterization studies, choose physical options) 

► Simulation: running the model, either in history matching simulations or in 
forecast simulations 
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Reservoir simulation models 



► What is Reservoir Simulation? 

► What is a Numerical Simulator? 

► Reservoir Simulation 

• Is an approach that enables the reservoir engineer to analyse phenomena that could not 
be described by simple "hand calculations" or "traditional calculations", due to their 
complexity 

► A Numerical Simulator 

• Is a numerical model translating the image of our understanding of the reservoir into 
mathematical equations 
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Reservoir simulation models 


► A reservoir model is a mathematical representation or translation of the 
image that geoscientists and reservoir engineers have of a reservoir 

► Numerical simulations allow reservoir engineers to understand processes that 
are in action in the reservoirs and to predict reservoir behaviour under 
various production scenarios 

► Human factor is therefore an important aspect of reservoir simulation: 

• A potential risk would be that users could consider their results as "representing the 
reality"... 

• Users experience is therefore an essential aspect in reservoir simulation 
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Reservoir simulation models 



► A model is a simplified representation of the reality which allows the study of 
physical phenomena 

► A model may be viewed as a transfer box receiving inputs and delivering 
outputs but: 

• Some principles must be clearly recognized: 

• The model must not be used as a black box 

• The output data cannot be better than the input data 

► To be useful in reservoir studies, a model must be workable (size, time and 
cost) and trustable (acceptable accuracy) 
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Types of simulators 



► Mass balance simulator (ex.: MBAL) 

► Numerical simulator (ex.: ATHOS, ECLIPSE, VIP...) 


• Single phase 

• Black Oil 

• Compositional 

• Chemical 

• Thermal 

• Fractured 
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Reservoir Model Selection 



► Not every question needs in-depth modelling 

► The model selection will depend on: 




Scope 
Process 
Available data 

Resources (time, budget, software) 



TIME 

AVAILABLE DATA 
COMPLEXITY 
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Mass balance simulator 



OHIP, (j), K, NTG, Swi 



► Very simple model 

► Material balance equation 

► Checking data consistency 
► Determine fluids in place 

► Estimate aquifer type and size 
► Production prediction 



THE RESERVOIR IS REPRESENTED BY A TANK WITH 
A UNIQUE POROSITY, PRESSURE, SATURATION 
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Mass balance simulator 
Material balance equation 



N.BOi = (N - Np)Bo + [(N.RSi - Gp ) - (N - Np)Rs]Bg 
+771. N. Bo t f 0 - -lj + We-Wp.Bw 


Introducing Rp definition: 



Gp 

Np 


Np[Bo + ( Rp — Rs)Bg] = N[(Bo — Boi ) + (Rs t — Rs)Bg] 
+ 777 . N. Bo t ( 0 - -lj + We-Wp.Bw 
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Mass balance simulator 


► Material Balance is suggested as a necessary step prior to carrying out a 
simulation study. 

► Material balance will always enable the drive mechanisms to be identified. 

► A material balance study can provide OHIP and drive mechanisms as inputs to 
simulation. 
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Numerical simulators 
Types 

► There are several types of numerical simulators 

• Black oil: One porous medium + black oil functions 

• Compositional simulator: One porous medium + EOS 

• Dual porosity simulator: Two porous media (matrix and fractures) 

• Thermal simulator: Flow + Energy equations. 

► A numerical simulator also integrates specific features such as: 

• Grid geometry (ID, 2D, 3D, CPG...). 

• Time step management. 

• Numerical schemes (discretization of equations). 

• Resolution methods. 

• Pressure drops through tubing and surface network. 

• Well schedule. 
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Numerical simulators 
Basic equations 

► Reservoir (volumetric): dV p = C p V p dP 


► Fluids (volumetric): 


dp f = -C f p f dP 


► Fluids (compositions): 


X x i = Z 3 / , = l 


x i 


y.- 



► Darcy's law (one phase flow): 


A (P - p.g.z) = - 


Qp 

K.A 


Ax 


> Multi phase flow: 


HPi-Pi-g-z) = 



K.kr t .A 


Numerical simulators 
Basic parameters 

► Reservoir (volumetric) 

• gross volume: Vt 

• net volume: Vu 

• porous volume & pore compressibility: Vp, Cp 

► Fluid 

• saturation: So, Sg, Sw 

• density: ro, rg, rw 

• viscosity: po, pg, pw 

• compressibility: Co, Cg, Cw 

• compositions: xi, yi, Ki 

► Flows 

• permeability: K 

• capillary pressure: Pcwo = Po - Pw , Pcgo = Pg - Po 

• relative permeability: krow, krw, krog, krg 


Numerical simulators 



► Space discretization 
Material balance equation per block 
► Properties distribution per block 



AX 

< > 



material balance in the block T 
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Numerical simulators 



material balance in the block T 


L 



time t 


( flow in) — ( flow out) ± ( production or injection term) = accumulation term 


l 



T Am-i 

Yu Qin ~ A 

^ Qout i Qi ~ 


MATERIAL BALANCE FOR EACH COMPONENT: 

BLACK OIL MODEL: OIL, WATER AND GAS 
COMPOSITIONAL MODEL: Cl, C2, C3, C4... Cn IFP Trainin g | 32 


Numerical simulators 
Basic equations: The flow term 



n K ab A ab kr ab 

Qo b = °-( Po - Po a - Po -9 -Az ah ) 


B r 


9o 


p os = stock tan k oil density 

p 0 = reservoir oil density 

P°,Po = oil pressure of cells 'a' and 'b' 

A z ab = z a - z b , depths of cells 'a' and 'b' 

g = gravity accelerati on 



ab = permeabili ty 
A ab = flows sec tion area 
dis tan ce 




transmissivity 


between cells 'a' and 'b' 


kr 0 = relative permeabili ty 
p 0 = vis cos ity 

B 0 = formation volume factor 


of oil between cells 'a' and 'b' 


J-l,k 


i-1, j, k 


U.M b 


a 


i+l, j, k 


i,j,k 


i, j, k-1 


i,j+l,k 
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Numerical simulators 

Basic equations: The accumulation term 



M = 


Am l 0 =V l 




S l 


P‘ 0 ^r + f AP '° 


V 


AP 


AP 


\ 


J 


APo + fi'pl AS l 


V 1 = volume 
(f)' = porosity 
p' 0 = oil density 
S l Q = oil saturation 


of cell '/' 


i, j, k+1 


i. j-1. k 


i-1, j, k 



i+l, j, k 


i, j+1, k 


i, j, k-1 
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Numerical simulators 
Production / injection terms 



► Constant well production rate. If the well is specified to have a constant well 
rate of Qi at reservoir conditions: 

► Constant well bottom hole pressure. For a well producing or injecting at a 
constant bottom hole pressure, Pbhi , the well rate is computed the following 

TD 

M, ' B, 


where Pli is the productivity or injectivity index of the well. 

► The PI may be specified externally, based on tests of the well. If the well is in 
the middle of the grid block, one may assume radial flow into the well, with 
block volume as the drainage volume: 

pl 2x-k r h 
ln (r,/r w ) 
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Numerical simulator 
Black oil model 



► Black-oil simulation means that the fluid properties depend only on pressure. 
They do not depend on the composition of that fluid 


• Simplified fluid description PVT - (Black Oil) 

• Reservoir rock properties / Petrophysics (K, 0 , Cp) 

• SCAL(Kr-Pc) 

• Boundary conditions: wells and aquifers 

• Initial reservoir condition 
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Numerical simulator 
Black oil model 

► Main assumptions 

• 3 components (in standard conditions) 

• Oil & gas compositions are assumed constant in time 

• Oil, water & gas density are assumed constant in time 

• 3 phases (in reservoir conditions) 

• Oil is a mixture of "oil" & "gas" components 

• Gas corresponds to "gas component" 

• Water corresponds to "water component" 

► 3 flow equations 

• Oil equation 

• correspond to "oil component" contained in oil phase 

• Gas equation 

• correspond to "gas component" contained in oil & gas phases 

• Water equation 

• correspond to "water component" contained in water phase 
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Numerical simulator 
Black oil model: densities 



► One phase gas: The gas must be single phase in the reservoir, which means 
that crossing of the dew point line is not permitted in order to avoid 
condensate fallout in the pores. 

_ Pgs _ constant 


► 


One phase water: Strictly speaking, it means that Pres > Psat of the water, has 
a density described by: Pws constant 


► One phase oil: Oil must be undersaturated, which means that the Pres > Pb. 
In the Black Oil fluid model: _ . _ D 

Pos + Pgs • K s 


Po = 


B, 


for undersaturated oil, Rs is constant, and the oil density may be written: 

constant 


Po = 


B r 


IFP Training 38 



Numerical simulators 

Black oil: The gas component flow equation 



P 'g,s 


ab a ab 


K aD .A 


- ab 


kr ab kr ' 

g sr>b r»a , „ „ \„ab\ , K ' 0 / r>b 


, ab 


B rV s 


( p; - p; + .g .Az ' 10 ) + Rs -2- ( p> -P° + p a .g .Az“ ) 

Vo 


= Am« 


J l_ 


FREE GAS 


DISSOLVED GAS 


Am” = V a A{f{p° g .S a g +p gs .Rs.S a 0 )) 


Am 0 “ =y°(p 0 ”.S: +p 9 , s .Rs.S:)A^ + 


rf 


9 9 


aV; _ APS 


A 


V 


AP 


AP 


+ s o Ap ; + + 


j 
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Numerical simulator 
Unknowns 



► Let's consider a mixture of n components in a cell (3 components in black oil 
model) 

► Unknowns linked to multiphase flow (for each phase): 

• Pressure 

• Saturation 

► Unknowns linked to compositions and phase equilibrium are (for each 
component): 

• Number of moles in mixture 

• Number of moles in vapor 

• Number of moles in liquid 
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Numerical simulator 
Unknowns 


SOLVING THE MATERIAL BALANCE FOR EACH COMPONENT: 


OIL, WATER AND GAS 


V ^ 

r 1 /3 m, 

( 


^ Qout i Qi ~ ^ 

+ 1 


C _L C _L C — 1 
J o ' g ' * J w A 

P — P = Pc 

r o r w r 

P — P = Pc 

1 g 1 o 1 ^ go 


WE OBTAIN: 



Numerical simulators 
Reservoir simulation model 





Q0SF3S 


0 23S07 


0*83713 
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Numerical simulators 
Compositional model 

► Black oil assumption is not valid ! 



► n Componentes: Cl, C2, C3, ..., Cn 

► 3 + 3n Unknowns: Mixture mass, L & V fractions, compositions (xn, yn, zn) 

► 3 + 2n Fundamental relations: 

• Liquid and vapour fractions: L + V= 1 

• Mixture compositions: L-x n +V-y n =z n 

• Liquid and vapour compositions: 2>n=l & 

• Equilibrium constants: y n /x n =k n {P,T) 

¥ We will have as many flow equations as components since we have one 
continuity equation per component 
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Other reservoir simulators 



► Compositional 

• Miscible gas injection 

• C02 injection 

• Gas condensate reservoir 

• Volatile oil reservoir 

► N components 

► Gas/oil mass exchange (Equation of state, constants) 

¥ Chemical 

• Polymers 

• Surfactants 

• Alkaline injection 

¥ Multicomponents in oil and water phase 
¥ Fluid/Fluid and fluid/rock interactions 
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Other reservoir simulators 


► Thermal 

• Hot fluid 

• Vapor 

• Combustion 

► Black oil or compositional thermodynamics 

► Reaction kinetics 

► Energy equation (temperature) 

► Fracture 

• Double porosity / double permeability 


Numerical formulation 
Partial differential equation 


Radial flow, ID (well test equation) 


transient flow: 


1 8 ( dP 

r 


r dr 


V 


dr 


(f>- /u-c 


P = P : + 


q-ju 


An-K-h 


K 


■Ei 


dP 

dt 


r , 2 \ 

</>■ ju-c-r ' 

4K -t 


Ez(-x) = -| 


" e~ u 


du 


steady state flow: 


1 d ( dP A 
r — 

V dr 


r dr 
P = P„ + 


= 0 

(p.-pj 

In (r,/r„) 


ln(r/r„) 


Numerical formulation 
Boundary conditions 



► External boundaries: 

• No flux at the boundaries 

• Aquifer limits 

► Wells: 

• Constant well bottom hole pressure (Dirichlet BC) 

• Constant well production rate (Neumann BC) 
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Well's boundary conditions 



► Rate specified in a well in block i. After computing the pressures, the actual 
bottom hole pressure may be computed from the well equation: 

Qo,i = PI i( P o,i- P bh,i) / (Mo,i' B o,i) 

► Production wells are normally constrained by a minimum bottom hole 
pressure, for lifting purposes in the well. If this is reached, the well should be 
converted to a constant bottom hole pressure well. 


► 


Constant well bottom hole pressure. Since rate terms contain unknown block 
pressures. At the end of each time step, actual rates are computed by these 
equations: 


R o,i = 


PI: 


<Po,i B bh,i ) Q w,i ~ 


PI, 


A - Ax i(Po,r B o,i) 

and water cut is computed by: fwsj 


A-Ax i (M w , i -B wJ ) 

R\,i 


( B w,i B bh,i ) 


R'w,i +C l'oJ 
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Numerical formulation 
Difference equations: ID 




a 2 p _( </> -^-c ^ ap 
8x 2 ~ t K J 


EXPLICIT FORMULATION: 


IMPLICIT FORMULATION: 


Ax 


V 


K 


pu-2p;+pu pr At -p- 


j 


At 


pt+At _ ry pr+Ar . p 

-M+l fAi ^i-1 

Ax 2 


t+A1 ' " ,+At 

^ • — 

l K 


At 


!FP Training 


51 


Numerical formulation 
Stability of formulations 



► Using the von Neumann stability analysis, the stability requirement for the 
explicit formulation is: 



► Therefore, the time step size is limited by both grid block size and properties 
of rock and fluid. 


► For the implicit formulation, the von Neumann analysis results in: 


At < oo 


► Showing that it is unconditionally stable for all time step sizes. 
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Methods of solution 



► Explicit: impractical, 

► Fully implicit (FIM): most robust, but expensive, 

► Partially implicit (IMPES): it is less dispersive and sometimes faster than the 
fully implicit method. 


► Adaptive implicit (AIM) is generally the optimum approach: is a compromise 
between the fully implicit and IMPES procedures. 


In ECLIPSE simulator , the methods IMPES and AIM are used 
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IMPES 



► Assumptions: Transmissivities, coefficients and capillary pressures are 
evaluated at time t 

TX o,i + l,2( P o,i + l ~ P„,) + - P o,i) - R'o ,i = C pooA P o,i ~ K, ) + C LA S w,i - S w,i) 

PX w,i+ 1/2 [(Po,i+l _ P oA ~ i PC ow,i+l ~ Pc ow, / ) ]"*" PX w,i-l/2 [(^ 0 , 1-1 — P o, i) — ( Pc ow ,i-l ~ P C ow,i)] _ Q w,i 
- r 1 f rp _p t \ + r t ( Q _ C ( 

^ pow , z V o,i o,i ) ^ sww , z V ^ w , i ^ w , z / 

► Water and oil equations are combined in order to eliminate water saturation 
(Sw) as an unknown: 

a i- p 0 ,i-i + b i- p 0 , l +c l -P 0 , M = d l 

where a, b, c are fonctions of T,C and Pcow, d includes all the terms at time t 

► Using either the oil or the water equation, we solve explicitly for water 
saturation at t + Dt 
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Convergence Problem 



► Bad convergence of Newton iterations 

• too large Dt 

• convergence criteria too large (not enough iterations) 

• not accurate solution of linear system 

• physical reasons 

► If many step cuts 

• change from IMPES to FIM 

• decrease the Dt increasing factor 

► Physical reasons 

• cells with small PV 

• too large PI 

• too rapid change of rate 

• too small rate instead of Q=0 

• discontinuity in kr curves 

• very high Kv (fast segregation) 

• values outside the range in PVT tables 
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Remarks 



► Many assumptions are necessary to build a reservoir model so calculations 
are only approximate and a good engineering judgment is required to 
evaluate input data. 

► The solution of the flow equations consumes a large part of the computing 
effort in a simulation. As a consequence, selecting an effective equation- 
solving method for a particular problem is an important step in managing the 
cost and controlling the difficulty of a study. 
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Space and time discretization 


► Space discretization: 

• Grid definition 

• One value per grid cell for each data 

• Calculation of flow transfers between cells 

► Time discretization: 

• Time step definition 

• Calculation of production data per time step 

• Calculation of reservoir data at the end of time step taking into account constant 
production data and constant reservoir data during the time step 

• Data at the beginning of the time step = explicit schema 

• Data at the end of the time step = implicit schema 

► Simulation costs increase more rapidly than the number of cells as 

• Calculation per time steps are proportional to the number of cells 

• Time step decreases with the cell dimensions 
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Geometry 



Block -Centered Geometry 

O 2000 4000 0000 3000 



Corner 'Point Geometry 

0 2000 4000 GOOD 8000 



Gordon Adamson et al., 

Oilfield Review, SLB. Summer 1996 


Block centered geometry features 
flat-topped rectangular blocks that 
match the mathematical models 
behind the simulator. 

Corner point geometry modifies 
the rectilinear grid so that it 
conforms to important reservoir 
boundaries. 

Three dimensional grids are 
constructed from a 2D grid by 
laying it on the top surface of the 
reservoir and projecting the grid 
vertically or along fault planes 
onto lower layers. 
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Spatial discretization: GRIPPING 
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Time discretization 



► Constraints: 


AT. < A T < AT 

min max 

A P - max( A P t ) < A P K 
AS - max( AS ( ) < A S n 


► Calculations: 

At „ + i = a - At „ 

being a > 1 

if 

AP „ + 1 > A ^max 

■ AP ntl / AP max 

if 

A ^n+1 > A ^max 

A tn S + i = A t„' A S„ + ,/AS^ 


At » + i = min ( A Ci. At „ s + i) 


if At » + 1 > At m,x 

if At » + 1 < At ml0 


A t„ + 1 = A *m,x 

A U = A t ni „ 
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Time discretization 


► Production parameters are assumed constant per time step 

• Flow rates 

• Oil, gas , water, aquifer influx 

• Pressures 

• Oil, gas, water 

• Bottom hole flowing pressure 

• Well head pressure 

• Saturations 

► Flow simulator calculates for each time step and each grid block 

• For each phase 

• 1 pressure 

• 1 saturation 

• n concentrations 

• For each component 

• 0 to 6 exchanged flow with neighbouring cells 

• 1 injection or production flow for cells connected to a well or to an aquifer 
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Simulator work flow 
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Generalities 


CXu 

\ 

► Basic equations 

• Continuity equation 

• Darcy's law 

• Equation of state: it will depend on the model used 

► Models 

• Black oil: Produced oil and gas in standard conditions have constant compositions 

• Compositional model: use of an EOS to describe the fluid (p) 

► Space discretization 

• One value per grid cell for each data (p, K, S, P) 

• Calculation of flow transfers between cells 

► Time discretization 

• Production data are calculated at each time step 
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The need for reservoir simulation 



► The need for Reservoir Model comes from the fact that: 

• It is desirable to "see, or to predict the future" of a reservoir or a field, under various 
production mechanisms before deciding on a development plan and investment strategy 

• It is necessary to develop techniques to simulate (assume the appearance of) reservoir 
behavior under likely production schemes 

• It is desirable to integrate all the available reservoir data into a model 
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Reservoir simulation objectives 


► In other words, reservoir simulation objectives are: 

• To establish reliable production forecasts 

• To help decision making on development plan, appraisal drilling, infill drilling, 
workovers... 

• To improve reservoir description 
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Incentives for running a flow simulation 


► Main drivers: 

• To overcome simplifying assumptions (used in hand calculations) 

• To consider realistic problems (development plan, production history...) 

• To take into account all the available data and to better understand all the interactions 
between them for a given scenario 

• To perform sensitivity to unknown parameters 

► Main use: 

• Help to decision, in particular at the end of appraisal phase and during pre-development 
screening 

• Optimise the production profile and the final recovery after production start-up by 
improving the reservoir description 

• Establish reliable production forecasts 
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Incentives for running a flow simulation 



At all stages of a field development, dynamic reservoir simulation may help 

► Non producing reservoir 

• To identify best production strategy (depletion, injection...) 

• To define the well number and their architecture 

• To design surface facilities 

• To evaluate risk and profitability of the development project 

• To optimize CAPEX 

► Appraisal 

• To identify key dynamic uncertainties 

• To define appraisal needs (contacts, faults, facies variations) 

► Producing reservoir 

• To integrate all history production data 

• To improve reservoir description 

• To optimize oil production and recovery (infill, W.O, IOR...) 

► At all stages of the development 

• To establish reliable production forecasts 

• To evaluate (remaining) reserves 
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Time 
(not to 
scale) 
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Running a flow simulation 
Development project 

► Some examples of needs 

• To define the best production strategy (natural drive, water injection, gas injection) 

• To define the number, type and architecture of the wells 

• To establish production profiles (oil, water, gas) that will be used by economists and 
architects 

• To evaluate the risks associated to the project 
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Running a flow simulation 
Development project 

► Study approach 
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In this case, static and dynamic uncertainties are high 

One of the main objectives will be to evaluate how the different development scenarios 
will impact the production profile 

The need will specifically be to analyse the influence of certain reservoir parameters 
(permeability, fluid compressibilities, mobility ratios, heterogeneities...) 

This will necessitate a great number of prediction runs 




Running a flow simulation 
Producing reservoir 

► Some examples 

• To establish reliable production forecasts: Rates, WCut, GOR, Pressures 

• To optimise oil production and recovery whilst respecting all production constraints 
(pressures, surface facilities capacity, contract...) 

• To localise remaining oil, by-passed oil 

• To evaluate the need for infill drilling 

• To evaluate Work-over programs to optimise oil production or reduce water production 
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Running a flow simulation 
Producing reservoir 

► Approach 

• The History Match helps to reduce the static and dynamic uncertainties 

• The most influent parameters must be identified and adjustments to the reservoir 
description will then improve the reservoir behaviour description 

• That History Match process will necessitate a great number of simulation runs 
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Key points to keep in mind 



CXu 

\ 

Why run a reservoir simulation model 

► To define appraisal needs (contacts, faults, facies variations) 

► To evaluate how the different development scenarios will impact the 
production profile 

► To optimise oil production and recovery whilst respecting all production 
constraints (pressures, surface facilities capacity, contract...) 

► To identify key dynamic uncertainties 

► To assess the risks associated to the project 


Help to decision, in particular at the end of appraisal phase and during pre- 
development screening 
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General workflow 
and data review 


Reservoir Simulation Planning 



► Reservoir simulation study duration: from weeks to years 

► Necessity to plan carefully the study to give correct results in time, before to 
take decisions for the field management: 

• Problem definition 

• Data review 

• Data acquisition 

• Approach selection 

• Reservoir characterization - Build geological static model 

• Upscaling to generate dynamic reservoir simulation model 

• Computing support 

• Initialisation 

• History matching 

• Prediction 

• Reporting 
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Problem Definition 


► Determine present reservoir performance and associated operating issues in 
order to define future performance 

• Collect information 

• Identify the problems 

• Determine the objectives 
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Data review 



► Data generally must be reviewed and reorganized because: 

• They come from several sources of information 

• The review evidences inconsistencies and gaps 

► The objective of data review is to: 

• Decide if the data quality is sufficient to conduct the study 

• Plan data acquisition for the study purpose 
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Data acquisition 
Main Types of Data 



TYPE 

Geometry 

Petrophysics 

Description 

Shape 

Porosity 


Thickness 

Permeability 


Dimensions 

Compress. 


Faults 

Pc's 


Contacts 

Kr's 



Swi Sorw 

Origin 

Seismics 

Logs, cores 


Geology 

Laboratory 


Well tests 

Tests 


Fluids 

Wells 

Production 

Bo, po... 

Coordinates 

Qo = f(t) 

PS - ” 

Completion 

Qg = f(t) 

Bw, pw... 

PI 

Qw = f(t) 

Separat. 

Cond. 

Laboratory 

Well report 

Field 

Abacus 

Tests 
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Study approach 



► What is the most suitable model to answer the problems: 

• Model dimensions (ID, 2D, 3D) 

• PVT type (Black oil, compositional) 

• Special physical options (double media, thermal, chemical...) 

► Factors influencing the study approach are: 

• Availability of simulators and computing environment 

• Time, money and manpower 
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Model design 

► Model design is influenced by: 

• Type of process to model 

• Complexity in fluids mechanisms 

• Objectives of the study 

• Quality of data 

• Needed level of details 

• Budget 
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Grid Parameters 





REALITY 


GEOLOGICAL MODEL 
(OR GEOMODEL) 

for static simulation 


Phi 

K 

Sw 

^Pc 

Hu/Ht 

- A quantitative synthesis 

- One model for each parameter 

- Used to: 

-Understand spatial hydrocarbon 
distribution 

-Accumulation computation 


RESERVOIR MODEL 
for dynamic simulation 


Not to predict the reservoir 

But to anticipate 
the dynamic behavior 


History match 


► In general two phases are distinguished: 

► Initialisation 

• Model data checking 

• Calculation of fluids in place 

► History match 

• Run with imposed rate 

• Calculation of P, Fw, GOR 

• Comparison with field measurements 

• Adjustment of parameters 

• Modification of input data 

• New run 


History match 


► Objective: 

• Reproduce with the model the measured evolutions of pressure, BSW and GOR by well, 
by zone or for the entire field 

► Difficulties: 

• Uncertainties on fault and flow barriers network; parameters are dependent 

► Rule of thumb: 

• Predictions are reliable on a period twice the production period 
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History match - Principle scheme 
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Predictions 



► Predictions can be run once history match has been completed. The objective 
is to predict the future performance of the reservoir under different 
exploitation scenarios 

► Results of major interest are: 

• Future oil performance 

• Future water-cut and GOR evolution 

• Pressure evolution 

• Fluid contact evolution 

• Work-over requirements 

• Drilling requirements 

• Surface facilities requirements 

• Estimate the recovery factor 

► The analysis of results must be done by comparison of the different cases 
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Reporting 


► Give always a clear and concise report with: 

• Statement of the objectives of the study 

• Description of the model 

• Presentation of the results 

• Conclusions and recommendations 

► In most of reservoir studies give both: 

• Managerial report (read by many people) 

• Detailed report (read by only a few people, but essential for future update) 

► Don't forget to save data, results and reports. 
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Key points to keep in mind 



<Xu 

"\ 

► Simulation models are only tools 

► An integrated reservoir study requires a multidisciplinary approach that integrates 
complementary techniques 

► Use material balance and simple tools to check your data 

► Watch out during data analysis and input 

► Check the geologic validity of the data 

► Check result validity 

► A good model is one that provides an accurate flow representation, with: 

• An adapted 3D grid 

• A realistic set of petrophysical parameters 

• A successful upscaling 

► The model is reliable at time t, but at t + At? 


Be very patient! 
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Reservoir simulation workshop 
Simulator presentation 



Simulator presentation 


► Introduction to ECLIPSE 93 

► ECLIPSE input data and reservoir simulation 99 

► ECLIPSE keywords in well definition 117 

► ECLIPSE output 

Visualization of results 127 
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► The ECLIPSE simulator suite consists of two separate simulators: 

• ECLIPSE 100: Black-oil simulator 

• ECLIPSE 300: Compositional simulator 

► Both programs are written in FORTRAN77 and operate on any computer with 
a FORTRAN77 compiler 

► ECLIPSE 100 is a fully-implicit, three phase, three dimensional, general 
purpose black-oil simulator with gas condensate options 

• Oil and Gas phases are represented by one component 

• Composition of oil and gas components are assumed to be constant with pressure and 
time 

• The properties of those components can however change with Pressure and 
Temperature 
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ECLIPSE 



► ECLIPSE 100: 

• 3-D, 3-phase fully implicit Black Oil 

• Advanced local grid refinement, Polymer flooding 

• Pipeline network and wellbore friction model 

• Gas field operations and solvent model 

► ECLIPSE 300: 

• N component compositional and black oil model 

• Cube Equation of State 

• K-value or Black Oil treatment for PVT 

• Full coupling with network and geochemical models 
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How does ECLIPSE work? 



► ECLIPSE is a batch program, running with a single input data file 

► This Data file contains a complete description of the model: 

• Structure, reservoir description, fluid and rock properties, initial conditions, wells with 
their flow rates, surface facilities and constraints 

► The input file is structured by a serie of KEYWORDS and comments 

► This file is typically divided into 8 sections 
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Section-header keywords 


► Eclipse Model: the input data file is an ASCII text file organized into 8 sections 

RUNSPEC 

•• -j 

EDIT 

PROPS 

REGIONS 

SOLUTION 

SUMMARY 

SCHEDULE 

► The sections must be specified in the shown order. 

► It is recommended that the body of sections which are not frequently changed be 
held in separate files which are included in the data using the INCLUDE keyword. 

P A data record has to be ended with a slash [/] 
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ECLIPSE input data 
and reservoir 
simulation 


Section-header keywords 


RUNSPEC 


Title, problem dimensions, switches, phases present, components etc. 


GRID 


The GRID section determines the basic geometry of the simulation grid and various rock 
properties (porosity, absolute permeability, net-to-gross ratios) in each grid cell. 


EDIT 


Statu 

is: 

OPTIONAL 


Modifications to calculated pore volumes, grid block centre depths and transmissivities. 


PROPS 


Status: 


Tables of properties of reservoir rock and fluids as functions of fluid pressures, saturations 
and compositions. Contains the equation of state description in compositional runs. 
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Section-header keywords 


REGIONS 


Stati 

js: 

OPTIONAL 


■ Splits computational grid into regions for calculation of: 

■ PVT properties (fluid densities and viscosities): PVTNUM 

■ saturation properties (relative permeability and capillary pressure): SATNUM 

■ initial conditions, (equilibrium pressures and saturations): EQLNUM 

■ fluids in place (fluid in place and inter-region flows): FIPNUM 

■ If this section is omitted, all grid blocks are put in region 1. 


SOLUTION 


Status: 


Specification of initial conditions in reservoir may be: 

■ calculated using specified fluid contact depths to give potential equilibrium 

■ read from a restart file set up by an earlier run 

■ specified by the user for every grid block (not recommended for general use) 

This section contains sufficient data to define the initial state (pressure, saturations, 
compositions) of every grid block in the reservoir. ipprraining I m 


Section-header keywords 



SUMMARY 



■ Specification of data to be written to the Summary file after each time step. Necessary if 
certain types of graphical output (for example water-cut as a function of time) are to be 
generated after the run has finished. If this section is omitted no Summary files are created. 



■ Specifies the operations to be simulated (production and injection controls and constraints) 
and the times at which output reports are required. Vertical flow performance curves and 
simulator tuning parameters may also be specified in the SCHEDULE section. 
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Keywords in sections 
RUNSPEC section 



TITLE 


DIMENS 


title 

number of blocks in X,Y,Z directions 


OIL, WATER, GAS, 
VAPOIL, DISGAS 


the active phases present 


FI ELD/METRI C/LAB 



unit convention 

well and group dimensions 

indicates that input files are unified 

indicates that output files are unified 
start date of the simulation 

data checking only, with no simulation 
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Keywords in sections 
GRID section 



TOPS 


depths of top faces of grid blocks for the current box; data is taken from 
Structure map, and geological model from PETREL 


DX, DY, DZ 


X,Y,Z-direction grid block sizes for the current box; data is taken from Isopac 
map, and geological model from PETREL 


PERMX, PERMY, PERMZ 


X,Y,Z-direction permeabilities for the current box; data is taken from Isopac 
map, and geological model from PETREL 


grid block porosities for the current box; data is taken from Isopac map, 
and geological model from PETREL 


Grid and properties distribution from PETREL model 
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Keywords in sections 
PROPS section 




water relative permeability and capillary pressure as 
functions of S w 

oil relative permeability as a function of S 0 
in three phase system 

gas relative permeability and capillary 
pressure as functions of S g 


saturation tables 
from special 
core analysis 


FVF and viscosity of live oil as functions of pressure and R s 
FVF and viscosity of wet gas as functions of pressure and R v 



reservoir fluid 
properties from 
PVT analysis 


FVF, compressibility and viscosity of water 
stock tank fluid densities 


rock compressibility 
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Keywords in sections 
PROPS section 




Saturation tables from SCAL 


PVT file might be generated in PVTi using a PVT report 
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Keywords in sections 
REGIONS section 




fluid-in-place regions 


saturation table for drainage process regions 
saturation table for imbibition process regions 


equilibration regions 
PVT data regions 


Rock type distribution 
from PETREL model 
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Keywords in sections 
SOLUTION section 



EQUIL 


fluid contact depths and other equilibration parameters 


RESTART 


name of the restart file 


RPTSOL 


report switches for SOLUTION data 


Model initialization: 

Data from RFT, logs and MBAL 
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Keywords in sections 
SUMMARY section 




Field quantities 
Group quantities 
Region quantities 
Well quantities 
Connection quantities 
Block quantities 
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Keywords in sections 
SUMMARY section 



FOPT 

FOPR 

FGOR 

FWCT 

FWIR 

FWIT 

FOIP 

FOE 

FPR 

FORFR 

FORFF 

FORFW 

FORFE 


WOPT 


WOPR 


WGOR 


WWCT 


WWlH 


WWIT 


FQRFG 

FORFS 


Field or Well Oil Production Total 
Field or Well Oil Production Rate 
Field or Well Gas-Oil Ratio 
Field or Well Water Cut 
Field or Well Water Injection Rate 
Field or Well Water Injection Total 
Field Oil in Place 
Field Oil Efficiency 
Field Pressure 

Field Fraction Total Oil Produced by Rock Compaction 

Field Fraction Total Oil Produced by Free Gas Influx 

Field Fraction Total Oil Produced by Water or Gas Influx 

Field Fraction Total Oil Produced by Oil Expansion or Gas Expansion 
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Keywords in sections 
SUMMARY section 



WBHP 


WBP5 


WPI5 


WOPRH 


WWCTH 


WMCTL 


TCPU 


WTHP 


WBP9 


WPI9 


WORTH 


WGORH 


Well Bottom Hole Pressure or Tubing Head Pressure 
Well 5 or 9-Point Pressure Average 
Well Productivity Index Based on Value of WBP5 or WBP9 
Well Oil Production Rate or Oil Production Total History 
Well Water Cut or GOR History 

Well Mode Control 
Time of CPU Usage 
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Keywords in sections 
SCHEDULE section 



TUNING 


WELSPECS 


COMPDAT 


WCONPROD 


WCONINJE 


WCONHIST 


DATE 


time step and convergence controls 

introduces a new well: name, wellhead position, bottom hole reference depth, etc 
specifies the position and properties of one or more well completions 
control data for production wells 

control data for injection wells 
observed rates for history matching wells 

advances simulator to specified report dates 


Wells and production data 
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Datafile structure 



(^Repeat counts^) 


■ In the data following a keyword, asterisks may be used to signify repeat counts. A data quantity 
can be repeated a required number of times by preceding it with the required number and an 
asterisk. 


■ There must be no intervening blank spaces next to the asterisk on either side. 


■ For example, the two following cases are equivalent: 



RPTSCHED 


00000000 
0 2 0 0 2 / 

0 0 

RPTSCHED 


11*0 2 2*0 2 

/ 


■ Any lines beginning with the two characters are treated as comments, and will be ignored by 
ECLIPSE. 


-- VARIATION 

OF INITIAL RS WITH DEPTH 

DEPTH 

RS 

RSVD 


8200 

1.270 

8500 

1.270 / 
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Datafile structure 




■ Certain items of data can be defaulted to a built-in default value. The keyword description will indicate 
when defaults can be applied. There are two ways of setting quantities to their default values: 

1. by ending a data record prematurely with a slash (/) the quantities remaining unspecified will be 
set to their default values. 

2. selected quantities positioned before the slash can be defaulted by entering n* where n is the 
number of consecutive quantities to be defaulted. There must be no blank space between the 
number and the asterisk. If there is only one item at a time to be defaulted, then 1* must be 
entered (an asterisk by itself is not sufficient). 


WCONPROD 

'PRODUCER' 'OPEN' 'ORAT' 20000 4* 1000 / 

/ 
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Running simulation: launch ECLIPSE 



Running simulation 



► While the model is running, status reports will be written to the file 
FILENAME. PRT. 


► When the run is finished, your folder contains the following files: 

✓ FILENAME. DBG 

✓ FILENAME. EGRID 

✓ FILENAME. PRT 

✓ FILENAME. INIT 

✓ FILENAME.SMSPEC 

✓ FILENAME.UNRST 

✓ FILENAME. UNSMRY 

✓ FILENAME.RSSPEC 

✓ FILENAME. INSPEC 
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Well definition 


WELSPECS : General specification data for wells 

WELSPECS 

--name group i j BHP_ref_dep phase 
‘PI’ 'PROD 1 20 7 2500 'OIL' / 

/ 

■ Well PI belongs to group PROD 

■ Well head is at i=20, j=7 

■ BHP reference depth of 2500 m. Defaults to depth of top-most connection 

■ OIL is the preferred phase (used only for PI output) 

■ Other items can usually be defaulted 
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Eclipse keywords 
Well completion 



COMPDAT 
--name i 

j 

kl 

k2 

status 


CF 

diam 

kh skin 

direction 

‘PI’ 

20 

7 

3 

3 

'OPEN' 


23.47 

0.15 

I 


‘PI’ 

20 

7 

4 

4 

'OPEN' 


6.14 

0.15 

/ 


‘PI’ 

20 

6 

4 

4 

'OPEN' 


8.25 

0.15 

/ 


‘PI’ 

20 

6 

5 

5 

'OPEN' 


94.70 

0.15 

520.3 2 

1* Z / 


/ 


■ PI is a deviated well crossing columns (20,7) and (20,6) completed in layers 3 to 5 

■ The CF have been calculated in SCHEDULE application and input in item 8 

■ the well bore diameter must be given 

■ kh, skin and direction of penetration may be given for information as in the last line above 
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Eclipse keywords 
Production Well Controls 



WCONPROD : Control data for production wells 

WCONPROD 

--name status control Qos Qws Qgs Qls Qfond BHP THP VFP 

‘PI’ 'OPEN 1 'ORAT 1 500 300 1000 1* 1* 50 15 II 

/ 


■ Well PI has a target oil rate of 500, subject to: 

■ Maximum water rate of 300. 

■ Maximum gas rate of 1000. 

■ Minimum BHP of 50. 

■ Minimum THP of 15 

■ VFP table 1 applies to this well. 
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Eclipse keywords 

Well Economic Constraints 



WECON : Economic limit data for production wells 

WECON 

--name Qomin Qgmin maxWCT maxGOR maxWGR action 
‘PI’ 50 1* 0.8 1* 1* ‘CON’ / 

/ 


■ Well PI is subject to: 

■ A minimum oil rate of 50 sbpd: connections of PI will be closed if this rate cannot be met. 

■ A maximum water cut of 0.8: the connections with the highest WCT will be closed 
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Eclipse keywords 
Injection Well Controls 



WCONINJE : Control data for injection wells 

WCONINJE 

-name phase status control Qs Qres BHPlim THPlim VFP 
‘II’ 'WATER' 'OPEN' THP 1000 1* 300 40 2 / 

/ 


■ Well II injects water at THP = 400 bar, subject to: 

■ Maximum rate of 1000 sm 3 /d. 

■ Maximum BHP of 300 bar. 

■ VFP table 2 applies to this well. 
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Eclipse keywords 
Other well controls 



WELOPEN : Opens and closes both wells and completions. 


WELOPEN 


name 

action 

‘PI’ 

‘OPEN 7 

'P2' 

'STOP' / 

‘P3’ 

'OPEN' 

‘P4’ 

'SHUT' 

‘P5’ 

'SHUT' 


opens PI, does not alter completions 
closes P2 but allows cross-flow 
0 0 0 / opens all completions in P3 
0 0 2 / closes all completions in P4 in layer 2 
0 0 0 4 6 / closes completion 4, 5 and 6 of P5 


/ 


■ Opening a well and its completions requires 2 lines. 

■ Possible well status: OPEN, SHUT, STOP or AUTO 

■ Possible completion status: OPEN, SHUT or AUTO 
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Eclipse keywords 
Other well controls 



WEFAC : Defines well activity or efficiency factor (downtime). 

WE FAC 

-name efficiency 

‘PI’ 0.95 / 

/ 


■ The well PI is 95% of the time on production, for example to model shut-down of 1 day 
each month. 

■ Well rates and pressures are calculated using full rate, group rates and all cumulatives 
account for the efficiency factor. 
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Eclipse keywords 

Timestep management and convergence 


TUNING : Sets simulator control parameters. 

Record 1: Time stepping controls 

Record 2: Time truncation and convergence controls 

Record 3: Control of Newton and linear iterations 


TUNING 

1 50 0.1 / Timestep = 1 day, upper limit of next timestep = 50 days, and 
minimum length of all timesteps 

/ 

1 * 1 * 25 / Minimun & maximum number of Newton iteration defaulted, 

and maximum number of linear iteration = 25 
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ECLIPSE output 
Visualization of 
results 




Running simulation: launch ECLIPSE 
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3D visualization 
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Dynamic reservoir simulation 
Data review: gridding 



Sommaire 


► Review en geological modeling workflow 141 

► Reservoir description - Gridding 169 

► Grid types 181 

► Main issues in reservoir description 207 
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Review en 
geological 
modeling 
workflow 


Field Integrated Study Workflow 
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Construction of the structural model 


► With the set of faults selected in the 
fault model 

► Within the model boundary 

► Pillar griding 

► QC 



IN DEPTH (m) 


► Zigzag faults 

► Interpreted surfaces 

• Top A 

• Top B 

• Top C 


► Other surfaces 

• Interpolation between interpreted 
surfaces, constrained to well tops 
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Depositional Environment A 
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Packages of sand in fluvial channels 


Shaly packages of sand in levees 

Shales packages corresponding to 
floodplains 


FLUVIAL SYSTEM 



Buchans Brook, 
http://rainbowfish.anqfaqld.orq.au 
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Depositional Environment B 
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Sequence definition 
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Sequence definition 
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Layering 


Reservoir geometry: stratigraphy 

► Zones here corresponds to stratigraphic units. The objective is to refine the 
previous layering by creating the sub-layering corresponding to the units that 
one want to study using well markers and the previous horizons. 

► The purpose is to define the grid on which the computations will be done. The 
averaging of properties is done when we upscale the well properties. Then, 
the grid is filled with average properties. 




zones 
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Layering 


► Goal: find an optimal number of layers 

► Exact placement of layers will have a large effect only in models with 
correlated heterogeneities similar in scale to the displacement scale. 

► Honouring geological features 

► Honouring dynamic observations (RFT, PLT) 

► Enabling to simulate the well completion 

• Vertical 

• Horizontal 

• Selective, commingled 

► Enabling reasonable simulation duration 
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MFS Well to well correlation 








Geological grid definition - Vertical layering 
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Geological grid definition - Vertical layering 
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Zone definition 


• Definition based on the stratigraphic model (MFS, TS, SB) 
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From conceptual to numerical model: 
Layering 



► The type of layering must account for the depositional setting 


• Parallel to a surface ? 

• Proportional ? 
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Rock Types definition (Fluid flow model) 





Rock type 

RT1 

RT2 

RT3 

RT4 

RTS 

Porosity 

> 15% 

12.5 to 15% 

10 to 12.5% 

8 to 10% 


Irreducible SW 

5% 

8% 

14% 

20% 

45% 

Kv/Kh 

0.8 

0.5 

0.1 

0.05 

0.05 



Curve Sw vs. Heigth above Free Water Level 
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► Zone by zone 

► With the layering adapted to the depositional setting 

• Choice of a simulation method for Facies 
• U. Rutbah : pixel-based 


Mul. Clast : object based 
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Static Reservoir Model 


vertical proportion curve 


variograms 
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Petrophysical modeling 
Simulation results 

► Definition of the Porosity distribution 

► Correlation law between Log(K) and Phi 





Hor permeability 
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Reservoir 
description - 
Gridding 




Reservoir simulation model 



1 Reservoir simulation model 


1 Numerical calculator 
+ 

1 Set of modelled data 
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Gridding 



► In a numerical model, space is discretized and the reservoir volume is 
described as a series of contiguous cells defined by a grid system. 

► The finite-difference formulation of the reservoir simulation problem requires 
that the reservoir be segmented into a pattern of mesh points or gridblocks 


► Important parameters are: 

• The grid type 

• The cell size 

• The number of cells 
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Reminder 


Numerical Simulator Basic Parameters 

► Reservoir (volumetric) 

Gross volume: Vt 
Net volume: Vu 
Porous volume: Vp 
Pore compressibility: Cp 

Fluids 

Pressure: Po, Pg, Pw 
Saturation: So, Sg, Sw 
Density: ro, rg, rw 
Viscosity: po, pg, pw 
Compressibility: Co, Cg, Cw 
Compositions: xi, yi 

Flows 

Permeability: Kx, Ky, Kz 

Capillary pressure: Pcwo = Po - Pw, Pcgo = Pg - Po 
Relative permeability: Krow, Krw, Krog, Krg 
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Reservoir description 
Required properties 



► For each cell in the model, the following properties are required: 

► Geometry; identification of log correlations 

• Cell dimensions (x,y,z) and depths 

► Properties 

• Porosity 

• Permeability 

• Net to Gross or net thickness (if not indicated, usually assumed to be equal to 1) 




Cell properties such as: 



PORO 

PERMX 

PERMY 

PERMZ 

NTG 



are averages defined at the centre of the cell 
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Reservoir description 
Geometry 

► Rectangular grids: defined aerially by x columns and y rows 




► All columns have the same number of rows 

► All the cells in one column have the same length Ax 

► All rows have the same number of columns 

► All the cells in one row have the same width Ay 


IF* 


► 
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Reservoir description 
Geometry 

► Depths are defined at the grid block centre 



Ax 


Ay 


o 


XY view 



► Cell volume is Vt = Ax . Ay . Az 
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Reservoir description 
Grid 




Grids I et II are equivalent 
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Consequences of gridding 


► Numerical dispersion 

► Drawback of numerical techniques that can cause errors in simulation of 
processes involving rapid changes in saturation 

► No satisfactory way to eliminate numerical dispersion 

► Using of a large number of gridblocks reduces numerical dispersion 
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Consequences of gridding 


► Grid orientation effect 


► P and Saturations are distorted by grid orientation 


► As a consequence calculated performance is affected by grid orientation 


► Effect is higher when displaced phase is less mobile than displacing phase 
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Gridding 


: grid orientation 



► Grid selection will result from a trade-off between 

► Good vertical and horizontal reservoir description for: 

• Structure 

• Faults 

• Heterogeneities: vertical - horizontal 

• Initial fluid interface: WOC - GOC 

► Accurate saturation front description 

• Near wells 

• Near fluid interfaces Water-Oil and Gas-Oil 

• Between injection and production 

► Account for well location: present and future 

► Optimization of the number of cells 
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Gridding rules 



► Choosing a 3D grid consists in finding a compromise between various 
constraints such as: 

• Geological constraints 

• Dynamic (flow, production) constraints 

• Simulation constraints 

► Choosing a 3D grid means deciding the following points: 

• Grid dimension (number of cells) 

• Grid orientation 

• Grid block size 
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Grid size selection 



► Wells are generally considered in the center of cells 

► Leave a few cells between wells (at least 2) 

► Avoid large size variations between cells: ratio 2 

► Limit the number of cells: 

• Matrix size f(number of cells)2 

• Linear system solution 30 - 50 % of CPU time 

• Time step decreases with cells size (stability criteria) 

► Grid system must be carefully selected at beginning of study 
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Gridding 

Constraints and factors 

► Increase of CPU time with active cells 

► Orientation 

► Geometry / Geologic boundaries: 

• Faults 

• Horizontal barriers 

• External limits 

► Avoid large size variations between cells: ratio 2 

► Cell sizes: shape ratio and progression ratio 
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Active and dead cells 



► Dead Cells 

• Their porous volume is 0 

• Off-reservoir cells, allowing to represent a reservoir with an arbitrary form within a 
Cartesian geometry (parallelepiped) 

• Vp = (Ax * Ay * Az * 0 * Hu/Ht) * MULPV 

► Active Cells 

• Their porous volume is > 0 

► Desactive cells with low Vp (<MinPoreV) 


!FP Training 


186 


Geological constraints 


► The grid orientation should be parallel to: 

• The major faults direction, so that the grid block limits follow the fault trace. It produces 
a better representation of the transmissivity reductions if faults are present, with partial 
or total sealing effect 

• The fracture direction, in order to represent the permeability anisotropy kx vs. ky 

• The direction of sediment supply, for a better representation of drain, or to align with 
the main axis of turbiditic supply 
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Choice of vertical discretization 



► Geological constraints: Good representation of the limits (layers, units), the 
barriers and the drains 

► Dynamic constraints: Good representation of the displacement mechanisms: 

• Fingering, gas tongue 

• Gravity drainage effect (segregation, gas injection) 

• Oil rim 

► Numerical constraints: fine discretization around wells: 

• Good representation of well completions 

• Coning 

• Horizontal well 

¥ Alternative solution: LGR (Local Grid Refinement) 
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Vertical discretization 

► If it appears necessary to refine in "Z" geological layers: 

• Anticipate it as much as possible with the geologists before grid construction 

• In the other case, work with the geologist, in order to introduce the corresponding 
heterogeneities 

• Check if the given kz data are intrinsic 
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Using LGR to model gas coning 



■PERMX 

4950 to 
4984 to 
50 1 8 to 



Kx=Ky=Kz=5000 mD 
Qo = 6000 m 3 /j 
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Space discretization: sugar box geometry 



3 main flow directions 


1 cell can communicate with 
6 neighbours 
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Spatial discretization: Gridding - Geometry 




Comer 'Point Geometry 



Gordon Adamson et al., 

Oilfield Review, SLB. Summer 1996 


► Block centered geometry features 
flat-topped rectangular blocks that 
match the mathematical models 
behind the simulator. 

► Corner point geometry modifies 
the rectilinear grid so that it 
conforms to important reservoir 
boundaries. 

► Three dimensional grids are 
constructed from a 2D grid by 
laying it on the top surface of the 
reservoir and projecting the grid 
vertically or along fault planes onto 
lower layers. 
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Block-centered grid 



► Ease of use, no need for a grid builder 

► Traditionnally link to finite difference scheme 

► Poor description of complex geometry 

► Mainly reserved to phenomenological model 

► Physical connection between cells must be specified 

► No difference between dip or fault (physical flow not always respected 
between numerically neighboring cells): 



fault 

1.1 

Oi 

1 ,2 Q f 2 

Q2 





z 



Ql = — ( P 2,l~ P l,l) 
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Q2=—(P 2 ,2-Pl, 2 ) 

Q 2 = O 
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Ss 


^v- 
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“t- 


195 


Corner point grid geometry 



► Require a grid builder 

► Difficult grid building step when converting a geological model to a reservoir 
grid model if a complex reservoir geometry 

► NNC => irregular matrix system when cells are not orthogonal => require 2 
pressures for flux computation otherwise numerical errors 



CPG grid 




Align all corner on coordinate lines 



!PP Training 


197 


CPG grid features 



Pinch Out 



Tilted fault 
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Example of CPG reservoir model 





OOOOOC 0.23500 0 47000 D 70500 0 34000 
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Example of CPG reservoir model 




F 


OdSn 



0.3*479 0 4095? 0 6143* O.SISlS 
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Avantages of CPG grids 



► More flexible (8 independent nodes): 

► In vertical direction: 

• Dip description 

• Pinch-out modelling and tilted faults 

► In horizontal direction: 

• Accurate Fault trace 

• One transmissivity multiplier per fault 

► Physical connections between cells with common surface area: pinch-out and 
large fault throw 
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Draw back of CPG grid 



► Needs a grid builder 

► Difficult grid building step when converting a geological model to a reservoir 
grid model in a complex reservoir geometry situation 

► NNC -> irregular matrice system when cells are not orthogonal -> require 2 
pressures for flux computation otherwise numerical errors 
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Complex connection scheme 



fault 



Z 
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Hybrid grid 


Horizontal well: orthogonal grid 


Vertical wells: radial circular grid 


Reservoir: CPG grid 


Transition grid PEBI (Voronoi or finite element) 



Reservoir description: main issues 



► Reservoir layering 

• Identification of log correlations 

• Identification of flow units 

► Population of the petrophysical model 

• Definition of Net to Gross, Porosity and Permeability per cell 

• Two main problems 

• Interpolate (or extrapolate) data in non recognised areas 

• Give an average for each parameter in each grid cell 

► Upscaling 

• Reservoir engineering data are not constant in one grid cell 

• Hence the necessity for each parameter to precise: 

• The parameter distribution in one the grid cell 

• How to replace this distribution by one average value 
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Reservoir layering: use of log correlation 




Top U. B . 


Top L. B 


Top A F Dolo 


MFS3’ 

MFS2’ 


MFS5 

MFS3 

MFS2 


MFS1 


!f=P Training 


209 


Reservoir layering: identification of flow units 




Pressure (psia) 



Reservoir layering: non adjacent layers 



► It is possible to have non adjacent layers 


Reservoir 1 (Layer 1) 
Inter reservoir (no layer) 
Reservoir 2 (Layer 2) 


► This possibility can be used when 

• There are interferences between the two layers 

• There is no oil displaced in the inter reservoir 

• Flows leaving one layer are identical to flows going into the other layer 

► Example: Permeable layers communicating through a water bearing low 
permeability layer 
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Reservoir description: net thickness and porosity 




Facies 1 (<1)1, Vul) 
Facies 2 (<|>2, Vu2) 
Non reservoir 


Cell Volume: Vt = Ax . Ay . Az 
Net Volume: Vu = Vul + Vu2 
Porous Volume: Vp = 4>1 . Vul + 4)2 . Vu2 

Net Thickness: Hu = Ht . Vu / Vt 
Porosity: 4 > = Vp / Vu 
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Reservoir description: averaging data 



AVERAGE SHALINESS = 10% 


























































































































































































































































































































































1 



























































!PP Training 


213 


Upscaling 



■ Geological model contains a very 
large number of grid blocks. This 
level of resolution is generally 
incompatible with the computing 
capabilities of numerical flow 
simulators. 


- To do so, fine grid blocks are 
grouped in aggregates: coarse 
grid blocks. The basic problem is 
to determine the equivalent 
properties of the coarse 
gridblocks. 



Fine grid 
(geological model) 




Coarse grid 
(reservoir model) 
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Upscaling: which property? 


► Additive variables: Swi, Sgc, Vsh, NTG... 

• Exact volume averages 

► Absolute permeabilities 

• Analytical means: often used, but approximate 

• Numerical methods to minimize the loss of information 

► Relative permeabilities, capillary pressures and PI 

• Not that easy, often disregarded 

• Numerical techniques 
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Permeability 



► The permeability is not a scalar, but rather a second order tensor. 

► As permeability is not an additive property, Keq cannot be derived from a 
simple analytical formulation such as the arithmetic mean. Flow rate or 
dissipated energy must be considered instead: 

• flow rates at the boundaries of the aggregate of fine gridblocks must be the same as 
those of the coarse gridblock. 

• dissipated energy must be identical for the aggregate of fine gridblocks and the 
corresponding coarse gridblock. 
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Algebraic methods: parallel flow 



PARALLEL grouping 
AP constant 

AP identical for the 2 layers 



K 

k. 


AP = Q 1 L 


A 


IKK 


= Q 2 L 


A = (Qj+Q 2 )L// 
l h 2 k 2 l (h 1 k 1 + h 2 k 2 ) 

Compute k e to get Q with AQ 1 +AQ 2 


AP=(Q,+Q 2 ) 7 L/ \ 

2j l(h 1+ h 2 )k e 


AP 


III 


k e = 


h 1 k 1 + h 2 k' 

h j + h 2 


ii 


ll 



2 

Qi 
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Algebraic methods: serial flow 



l Z 




AP, 


AP, 


r 


Compute k e to get Q with AP 1 +AP 2 


AP = 


Qju 


Ih 


r 


V 


K , A 
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K k 


2 J 


Q ju L x + Li 2 

Th k~ 


Q 


111 


SERIES grouping 
Q constant 

Flow rate conservation 


A P 1 = 


AP , = 


Q 

/ k x 
Q L 2J u 

Th~k7 


_ L 1 + L 2 
e L, U 


k i 


+ 


Harmonic average of permeability perpendicular to the flow direction 
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Algebraic methods: random isotropic flow 



• It is rigorous in the chess case, anisotropy: 


Q 



ln(K g ) = -'£ln(K i ) 
n 


Mean of Cardwell & Parsons: 


K cp =(K 


CP , min CP , max 


J/2 


It can be shown that the Cardwell & Parsons mean bracket the equivalent permeability 


K u <K 


CP , min 


< K < K 


CP ,max 


eq 


< K 


flow || planes flow - 1 - planes 
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Algebraic methods: Cardwell and Parsons 
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Harmonic averaging Arithmetic averaging 
by row 


Flow based method 



► Tensor: 


No flow _L to AP 

Qx 



























◄ 

v u 

AP X 

K= a/APx 

Kxy=0 


Diagonal permeability tensor 


Flow based method 


Linear or periodic boundary 
conditions 

4* 

ft ^ 






0 % 


















V U 

A 


Ap x Kt = Q/AP x 

K(y = Q/APx 


Complete permeability tensor 
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► Numerical estimators differ from analytical methods by a need to solve the 
diffusivity equation: 


div(kV-P) = 0 


► for each coarse block, solve in each direction the pressure equation using 
finite element or volume technique, new techniques are proposed in order to 
capture the influence of the neighboring heterogeneities: 
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How to choose the aggregation rate? 



► Try to minimize the loss of information 

► Do not include large drains or barriers in a coarse block: look at the facies 
proportion 

► Compare permeability histograms 

► Use single phase simulations 

• Pressure field depletion 

• Production curves comparison 

► Keep enough grid blocks between injectors and producers to preserve 
breakthrough times. 
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Modelling high perm streaks 



► Some reservoirs contain one to several high permeability streaks with order 
of magnitude higher permeability than average. 

► Reservoir surveillance can successfully identify such streaks in reservoirs most 
seriously affected. 

► Both static and dynamic models of affected reservoirs should reflect the 
existence of high permeability streaks. 

► Technological challenge is to represent each and every high permeability 
streak identified in reservoir 
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Key points to keep in mind 



CXu 

\ 

► In a numerical model, space is discretized and the reservoir volume is described as 
a series of contiguous cells defined by a grid system 

► The grid must take into account sedimentological and structural elements in order 
to explain reservoir anisotropy regarding fluid flow as well as dynamic information 

► The gridding integrates reservoir geometry and petrophysical properties 

► The grid orientation should be parallel to the major faults direction, so that the 
grid block limits follow the fault trace 

► Grid geometry: cartesian, radial, non-orthogonal 

► Grid types: Block centered. Corner point geometry (CPG), Hybrid 

► Main issues in reservoir description: 

• Layering 

• Properties distribution 

• Upscaling: aggregation rate and properties' averaging 
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Dynamic reservoir simulation 
Data review: PVT 



Sommaire 


► Fluid description 228 

► Black oil model 234 

► PVT tests 250 

► Compositional model 264 

► PVT MODELING exercise 268 
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Phase envelope of a mixture 



Oil reservoirs with 
dissolved gases 


Q) 


V) 

cy 


Pcc 


Pc 



Gas reservoirs without 
retrograde condensation 

O 


Critical point 


Cricondentherm 


Dry gas 


Zone O: No or poor contribution of 
dissolved gases 

Zone ©: Appreciable contribution of 
dissolved gases 

Zone ©: Retrograde with condensation 
of liquid in the reservoir 


Zone 0 : Dry or wet gas 


Tc Tcc 
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PVT data for numerical simulation 



► The PVT model describes: 

• fluids (oil, gas, water) 

• densities 

• compressibilities 

• viscosities 

• phase equilibrium (between oil and gas generally) 

► The PVT model is chosen according: 

• to the hydrocarbon type 

• heavy, medium oil 

• light oil, condensate retrograde gas 

• and the production mechanism 

• water injection, gas injection 

• miscible displacement 
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Hydrocarbons: surface & reservoir conditions 



P stock 
T stock 


P res 
T res 
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PVT data for numerical simulation 


The two main PVT models are: 

► BLACK-OIL (improved BO, Bi-components) 

• tabulated data 

• simplified PVT 

• for heavy - medium oil (Bo < 2), water injection, ... 

► COMPOSITIONAL 

• Equation of State (EOS) 

• more and more precise and CPU consuming according to the number of components 

• for light oil (Bo > 2), retrograde condensate gas, miscible displacement, N2 or C02 
injection,... 
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Fluid description: influence of reservoir temperature 



Black oil: Small variations of composition 



BLACK OIL MODEL: small variations of composition 
COMPOSITIONAL MODEL: recomended near of critical point 
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Black Oil Model 


► 3 phases: WATER, OIL, GAS 

► 3 components: WATER, OIL, GAS 

► Reservoir temperature is constant with time 

► OIL and GAS in stock tank conditions have 
constant compositions 

► OIL and GAS in reservoir conditions are 
mixtures of pseudo components OIL and GAS 

► Volumes in reservoir conditions depends on 
pressure and concentrations 


STOCK TANK CONDITIONS 


GAS 


GAS 


OIL 


OIL 


Pos Rs Pgs 


R S Pgs 
r os + Rs Pgs 


C oo = 


RESERVOIR CONDITIONS 

Pos 

Pos + Rs Pgs 
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Black oil model 







Hydrocarbons: surface & reservoir conditions 



► Formation Volume Factor 

Boi = Bob [l+Co(Pb - Pi)] for P >Pb 
Bob = Boi [l+Co(Pi - Pb)] for P >Pb 


Volume of gas 
Volume of liquid 


Volume of liquid in reservoir conditions 
Volume of liquid in stock tank conditions 

_ Volume of gas in reservoir conditions 

Bg = 

Volume of gas in stock tank conditions 


API = 


141,5 

d 


131,5 (d = oil density en g/cm3) 
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Main oil properties (Bo, Rs, Co) 



► Solution gas/oil ratio (Rs) 

V g std Volume of gas in standard conditions 

V std ref "" Volume of oil in standard conditions 

Rs quantifies the amount of gaseous components which are dissolved in the oil at 
reservoir conditions. 


► Compressibility (Co) 


• Co = 


1 

V 


< av 

dP 


)t 


1 (4§-°)t 


Bo dP 


• Co quantifies the changes of oil volume during depletion at reservoir temperature, 
above the bubble point pressure. 
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Black oil: compressibility calculations 



OIL 



Bg ARs 



Bo 

1 

Bo - ABo 


P P-AP 

AVo = -ABo + Bg ARs 

^ 1 ( A Bo „ ARs 

Co = Bq 

Bo { AP A P ) 


GAS 




Bg 


Bg + ABg 


P P-AP 

AVg = ABg 

Bg AP 
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Main gas properties 



Gas specific gravity 

► Gas specific gravity is the ratio of the gas density @ std. cond. to the air 
density @ std. cond. 

Gas y can be determined from the gas composition and Mw: 

• y = (gas Mw)/(air/Mw) with Mw = 28.996 g/gmol @ 60°F, 1 atm) 

► Definitions for gas: 

• Gas Specific Gravity: y (air = 1) = p gas /p air = Mw gas /Mw ajr 

y = Mw gas /28.978 

• Gas density: p gas = Mw gas /V air(molar) = Mw gas / 23.645 

Vaird mole of air) = 23.645 p gas = y* 1.225 kg/m 3 

• Air density: Mw air = 28.9784 p air =28.9784/23.645 = 1.225 kg/m 3 

!FP Training I 239 







Black oil model 

Typical behavior of oil properties 



•Rs 


•Bo composite 


•viscosity 



od 

o 

CO 
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Oil properties - Solution GOR 



Standard conditions 



200 Sm 3 gas 


0,8 Sm 3 oil 


Example: 

Rs = 200/0.8 = 250 Sm3/Sm3 


► Solution GOR (Rs): 

• Ratio between the volume of gas liberated at standard conditions and the volume of oil 
obtained at standard conditions 

Rs = Vg lib (Pstdjstd) / Voil (Pstd, Tstd) 

• Rs quantifies the amount of gaseous components which are dissolved in the reservoir 
oil. 

• Units: Sm 3 /Sm 3 , Scf/Stbbl 

• Typical values : 100 Sm 3 /Sm 3 (black oil), 300 Sm 3 /Sm 3 (light oil) 
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Oil properties - Production GOR 



Standard conditions 


Reservoir conditions 



200 Sm 3 gas 


0,8 Sm 3 oil 


100 Sm 3 
gas 


Example: 

Production GOR = 300/0.8 = 375 
Sm3/Sm3 


► Production GOR: 

• Ratio between the volume of gas produced (liberated gas + free gas) at standard 
conditions and the volume of oil obtained at standard conditions 

GOR = (Vg lib (Pstd, Tstd) + Vg free (Pstd, Tstd)) / Voil (Pstd, Tstd) 

• Units: Sm 3 /Sm 3 , Scf/Stbbl 

• Production GOR > Solution GOR 
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Oil properties - FVF 


Standard conditions 



Reservoir conditions 



200 Sm 3 gas 


0,8 Sm 3 oil 


Example: 

Bo = 1.0/0.8 = 1.25 m/m 


► Formation Volume Factor (FVF or Bo): 




Volume of oil at reservoir conditions (P,T) 
at standard conditions 
Bo (P,T) = Voil (P, T) / Voil (Pstd, Tstd) 
Units: v/v, m 3 /Sm 3 , bbl/Stbbl 


which must be produced to obtain 1 m 3 of oil 
Do Pg XR s+PosT 


Typical values: 1.25 v/v (black oil), 2.0 v/v (light oil) 


Po(pT) 
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Fluid description: PVT Region 



► Constant Saturation Pressure 

• All samples collected at initial state have the same saturation pressure 

• Measurements of Bo(P), Rs(P), po(P) for one sample give values of Bo, Rs, po 
everywhere in the reservoir at anytime 

► Variable Saturation Pressure 

• All samples collected at initial state have a saturation pressure varying with an 
understandable law (e.g.: saturation pressure depends on depth) 

• Measurements of Bo(P,Pb), Rs(P,Pb), po(P,Pb) for one sample give values of Bo(Pb), 
Rs(Pb), po(Pb) everywhere in the reservoir at anytime 


PVT region: constant saturation pressure 
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Initial state 
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Depleted reservoir 


PVT region: constant saturation pressure 
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PVT region: variable saturation pressure 





Initial state 


Depleted reservoir 
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PVT region: variable saturation pressure 


Ech2 
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Oil PVT study - Differential liberation 


► Determine: 

• Bubble point pressure at reservoir temperature (Pb) 

• Amount of solution gas as a function of pressure 

• Liberated gas properties: compressibility, composition, density, Bg 

• Oil properties as a function of pressure: density, Bod, Bobd, Rsd, 


reservoir temperature T ( 


T = 15°C 


Vt3 

Vt2 


Volume 


Vo2 

Vo3 


li \ 


3 bubble point pressure 





P2< PI 


P3 < P2 


P = 1 bar 


!FP Training 


251 


Oil PVT Study - Constant composition expansion 



► Determine: 

• Bubble point pressure at reservoir temperature (Pb) 

• Isothermal compressibility factor of oil at bubble point (Co) 

• Gas compressibility factor (Z) 

• Global volume of oil as a function of pressure 
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Oil PVT study - Separator test 



► Objective 

• Simulation of field test or process separation scheme 

► Main properties derived from this experiment 

• GOR 

• Process R si/ Process B oi 

• Gas composition at each step 

• Stock tank liquid composition and properties gas 



OIL /FPTraining I 253 


Pr - Reservoir pressure 

Pb = Bubble pressure 

Pwf = bottom hole flowing pressure 




'/////////A 

Fluid produced 


Composite PVT 



► Differential liberation corresponds to the conditions in the reservoir: 

• isotherm 

• progressive depletion 

• liberated gas removed 


► From the well to the stock tank, liberation conditions are different: 

• temperature varies (from reservoir to stock tank) 

• separation through 1, 2 or 3 stages separator chain 

• gas and oil in contact in the well, specific process through the separation chain 

► Two additional curves are needed: 

• volume contraction (Boc) 

• gas liberated (Rsc) 


► To describe the well + separators chain to tank process: this constitutes the 
composite PVT 
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Constant Composition Expansion {CCE) 



Pj > Psat 



P2 > Psat 



P3 « Psat P4 < Psat P5 < Psat 


STEP 1: CCE 




P3 < P2 


T ■ 15*C 



P » 1 bar 


STEP 2: DL 







Multistage Separation Test 


V stage 


^wf 

T r 

Bottom hole 
conditions 


GAS 




conditions 


' 2° stage 
flash 




OIL 


STEP 3: MULTISTAGE 
SEPARATOR TEST 







GAS 

Stock tank 
conditions 


OIL 


) 


- Fluid produced 



Bo and Rs composites, P > Psat 



► Composite formation volume factor. Bo 


Bo 


composite 


( P ) = Bo ( Psat ) 


Vo(P) 


CCE 


Vsat(P) 


CCE 


► Composite dissolved gas, Rs 


Rs 


composite / r)\ r> flash 


(P) = Rs T (Psat) 
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Bo and Rs composites, P < Psat 



► Composite formation volume factor. Bo 


Bo 


composite 


(P)=Bo DL (P)- 


Bo flash (Psat )_ 

Bo DL (Psat ) 


► Composite dissolved gas, Rs 


Rs composite ( P ) = Rs flash ( Psat ) - (j Rs DL ( Pi )-Rs DL (P j)- 


Bo ftash ( Psat ) 


Bo DL ( Psat ) 
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Composite PVT 
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...Which PVT to choose? 


► A fluid sample will follow a composite path between its original location in 
reservoir and its final destination at surface 

► First step 

• Fluid is moving in the reservoir above saturation pressure 

• Volume change versus pressure is identical to CCE 

► Second step 

• Fluid is moving in the reservoir below saturation pressure 

• Liquid composition versus pressure is identical to a CVD 

► Third step 

• Fluid has reached the well bore 

• Volume and composition surfaces are identical to a flash process 
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Compositional 

model 




Compositional model 



► n pseudo-components 

• Cl, C2C02, C3C6, .... 

► Characteristics of each pseudo-component (concentration, Tc, Pc, Vc, ...) 

► Volumes computed with an Equation of State 

► Phase equilibrium from 

• Equation Of State 

• Tables 

► Each pseudo component is present in reservoir gas and in reservoir oil 
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Main EOS in reservoir simulators 



► Redlich-Kwong 


► Soave-Redlich-Kwong 


p _ RT 

aT ~ l/2 

V-b 

V-(V + b) 


p_ RT 

o(T) 

V-b 

V-(V + b) 


► Peng-Robinson 


RT a(T) 

V-b V -(V + b) + b-(V -b) 


► Zudkevith-Joffe 
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► The PVT model is chosen according to the hydrocarbon type and the 
production mechanism 


► The two main PVT models are: 

• BLACK-OIL: 3 components OIL, GAS, WATER (simplified PVT with tabulated data) 

. COMPOSITIONAL: Equation of State (EOS) 

► PVT properties in BO model: 

• Formation volume factor Bo= f(P) and solution gas/oil ratio Rs = f(P) 

► PVT experiments in order to determine main fluid properties 

• Differential liberation: differential FVF and Rs, liberation GOR 

• Constant composition expansion: bubble pressure, compressibility 

• Separator test: composition, viscosity, density, FVF and GOR 

► Composite PVT is used to describe the process from reservoir to surface, 
the three PVT experiments are combined: DL, CCE and separator test 

V 
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We define the different PVT analyses (CCE, DL, etc...) and observations (liquid 
density, GOR, etc...) in the main data tree. 

The data-set and project files are designated as *.PVI 

The export files to the ECLIPSE simulator are designated as *.PVO 

The print files with the simulation results are designated as *.PVP 


The main fluid is always called Zl. We may choose the name of the other fluids 
of the system, such as those corresponding to mud filtrate, swelling fluid, 
produced by mixtures, after mud decontamination, etc... 
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ECLIPSE PVTi PVT analyses 

Many PVT analyses can be modelled: 



► Single-point experiments: 

• Flash 

• Multiphase flash 

• Bubblepoint 

• Dewpoint 

• Saturation pressure 

• Saturation temperature 

• Critical point 

► Pressure depletion 

• CCE 

• CVD 

• DL 


• Injection study 

• Swell (by mole fraction) 

• Swell (by GOR) 

• Vaporization 

• First contact miscibility 

• Multiple-contact miscibility 

• Multi-contact test 

• Separators 

• Composition with depth 


Any of these can be used, either to enter laboratory measurements and 
perform a regression to better characterize the fluid model, or as a prediction 
for a given fluid model. 
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Typical PVT modeling workflow 



1 . Analyse PVT report. 

2 . Enter fluid composition. 

3. Enter carefully selected PVT analyses from a PVT report. 

4. Plot basic fluid characteristics (fingerprint plot, phase envelop, etc...). 

5. Enter all existing PVT measurement results for the corresponding fluid. 

6. Run a simulation and plot observed vs. calculated results. 

7. Run regression for the improvement of the fluid model. 

8. Split the heavy pseudo-components 

9. Run regression for the improvement of the fluid model. 

10 . Group components when necessary, 
n. Run regression. 

12 . Export fluid properties. 
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PVTi: Getting started 



► Start the ECLIPSE program Launcher. 

► Click on the PVTi button. 

► Select the version and working directory as required. 

► Start the software by clicking Run. 



The main window is compartmentalized in various areas: 



(T) 

J yelp 


m PVTI 

file Edit yiew gun ytiUtie-s £2ptic 

a* y II EE ft )Kt ^ 5ft The menu area allows the user to use the tools of the simulation and take various actions. 
# # Gk ^ ^ ^ @ S i 9 SC O S 




* UNTITLED 




CD 


The data tree area provides 
a view of the different fluid 
samples of the project and 
the various experiments and 
observations defined for 
each sample. 






ub-plots area stores other plots 


that are not presently displayed and can be accessed 
and exchanged with the main plot by double-clicking. 


Navigation 


CD 


* 


The main plot display, where various graphs can be 
viewed and worked on. 


Legend 


This area presents the 
legend and caption as 
well as the main 
information about the 
simulation. 


Caption 



PVTi License Expiry Date 31^7^2015 
Project UNTITLED PVI initialised 
Reading RUNSFEC^HODSPEC 
Error exit £ro* RUN5PEOSY5TEM input 
No components del ined 


The log window displays the actions 
taken during the project as well as 
main results. 



m 


Feng-Robrison Jhree-parameter) 


INS 





Opening a new file in PVTi 



For this exercise, we will work with the report: ALWYN WELL 3/9a-N3, in order to build a fluid 
model for ECLIPSE Black oil simulation. 

• When ECLIPSE PVTi is opened, a window prompts the user to select the data-set file. Select "Annuler". 
You should get the window as shown in the figure when a new file has opened. 


SI PVTi 

E‘l« Ed ft flun ytHrtin GpHont Window Melp 

ilm & t 

la & a CJ ra 


= T 'a S3 
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ECLIPSE PVTi basic features and worfkflow for Alwyn 
3/9a-N3 PVT report 



► Two sets of separator oil and gas samples were analysed in the laboratory. 
One set was recombined using a corrected GOR to perform a complete PVT 
study (reservoir fluid), including: 

• P-V relation by CCE 

• single-stage separator 

• two-stage separator 

• DL 

• viscosity by CCE 


► The reservoir fluid, recombined using a corrected GOR of 145.2 Sm 3 /m 3 , has a 
saturation pressure of 257.5 bar(g), whereas the reservoir pressure is 445.0 
bar(g). There is no bottomhole sample for this well. 


!FP Training 


275 



Summary of ECLIPSE PVTi analyses and observations: 
Alwyn 3/9a-N3 PVT report 



The analyses and observations to be entered in ECLIPSE PVTi for this well, and 
the source and content of the tables in the PVT report are described. They will 
be entered as presented. Please go through the different analyses and refer to 
the PVT report to build your fluid model. 


► P-V relation by CCE (CCE1) TABLE VIII 

• Relative volume 

• Liquid density 

► Single-stage separator (SEPS1) 

• Liquid molar fraction 

• Vapour molar fraction 

• GOR 

• Stock-tank oil density 


]- TABLE X 
]- TABLE IX 
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Summary of ECLIPSE PVTi analyses and observations: 
Alwyn 3/9a-N3 PVT report 

► Two-stage separator (SEPS2) 

• Liquid molar fraction 

• Vapour molar fraction 

• GOR 

• Stock-tank oil density 

► DL(DLl) 

• Oil relative volume 

• Liquid density 

• Z-factor 

• GOR 

• Vapour molar fraction 

► Viscosity by CCE (CCE2) 

• Liquid viscosity 

► Saturation pressure (PSAT1) 

• Saturation pressure TABLE I 


y TABLE XII 
y TABLE XI 

- TABLE XIII 

— > TABLES XIV, XV 
TABLE XVI 
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Enter a new fluid 



To enter the fluid composition, use the Table X from the PVT report in order to create Zl. In the case, the 
"Fundamentals" window is not open, to enter the fluid composition: right-click on "Zl" go to Edit Mole 
Fractions... enter the composition of the calculated reservoir fluid. 


You can observe on Table X, that you have a pseudo-component C11+, enter the mole weight and the 
specific gravity for the C11+. Click on Apply and check that the total is 100%. 


Close the window! 



i~o- ES 


Pi — B 
e — B 
— B 
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Plot the main fluid characteristic plots... 



The reservoir fluid, called by default Zl, appears now in the data tree. Look at the main characteristics of the fluids. 

a. Go to View / Sample / Fingerprint plot in the menu. You should observe the plot as in the figure, showing the 
composition vs. the mole weight of each component. Another alternative is to right click on Zl and select 
Fingerprint plot directly. 



Plot the main fluid characteristic plots... 



b. Do the same to observe the phase plot (View / Sample / Phase plot, or alternatively right click on Zl, or grab 
and drag Zl to the main graphic window). The plot should appear as shown in the figure. 



Note : Each time a new window opens, the previous one from the main plot area goes to the sub-plots area. You 
may navigate between these by selecting one and double-clicking. 
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Main windows for fluid characterization... 



Navigating through the main windows, you can see the main information and options and where they can be 
changed. Feel free to navigate through them! 

Units: To view and select the units, click on Utilities / Units. 

The units system we use for this exercise is Metric. Check the units! 



Mole fractions: To observe the fluid composition, go to Edit / Sample / Composition. The compositions of all the 
samples are displayed in this window. This is where you entered the fluid composition. To view compositions 
fluid by fluid, right-click on one of the fluids and select Edit Mole Fractions. 


y tifeta Fractions for Sample Zl 



28-1 



Main windows for fluid characterization... 



EOS: Go to Edit/ Fluid Model / Equation of State. 

You can select from four different EOSs and four different viscosity equations. 
The default selection is PR3 (Three-Parameters Peng-Robinson EOS) 
with correction term and Lohrentz-Bray-Clark for the viscosity correlation. 



♦ PR3 Peng-Rotwon 

SRK3 3-Pvamdef Soave-Re*cb-Kwoog 
SW: Sctwkdt Wwud 


PR ; 2-Parameter Peng- Ro bin so n 
SRK : 2- Parameter Soave ■ Rtdbch ■ KwOmg 
RK; Rtefrch-Kwoog 
Z J Zud^e-vtch-Joffe 


These are the best adapted in our case, don't change them! 


S Correction lam for PR, PR3 & SW EOS 
Viscorty Conatalion Type 
* Lrtireni-&ayOark 


P*dWMf> 

Aagt#s-P«t«n«n 
Lotveni-Biay-Qark (MocMiecf) 
Rwet EOS Parameters to d*f»A 


| OK | 


Components: Go to Edit/ Fluid model / Components and go through the different tabs. Names displays all the 
components and their types, whether they are defined by library, characterized or user defined. The type can be 
modified from this tab. In Complete the full characteristics of the each component are displayed (mole weight. Pc 
and Tc, O a and O b , etc...). This tab is for viewing purpose only. In Library are displayed the library characteristics 
of any component type defined as "Library" in the Names Tab. In the User tab, user defined characteristics can be 
entered. As there are no user-defined component types for this fluid model, it is empty. In the Characterization 
tab, the pseudo-components (here C11+) are defined as per chosen critical properties and acentric properties 
correlations. 

Check the different options but leave them as the default ones (Kesler-Lee) when you leave the window... 




Adding experiments... 


a. Go to Edit / Experiment.... A window opens where you can add a new experiment or edit the existing 
ones. Click on Add... 

b. Let's begin with a CCE Experiment: Once you have clicked on Add... then Pressure Depletion / Constant 
Composition Expansion... 

c. On the PVT report, you find The Constant Composition Experiment for the reservoir fluid "Zl" (Table VIII) 

d. In the General Tab you choose the fluid sample, Next... specify the fluid and temperature then Back... 

e. In the Observation Tab add the pressure data required for the experiment and select the column 
heading for the observations. Choose "Pressure", "Relative volume" and " Liquid density", then click on 
Next... fill in the table and click on Next... 

f. The window goes back to the General Tab. 

g. Add the experiment One stage separator (SEPS1): Edit / Experiments... /Add /Separators... then fill in 
the items: 

General Tab : Zl. 

Observation Tab: Define Temperature, Pressure, Liquid density, GOR (Table IX). 

Components Tab : Toggle 'Liquid molefrac.' and 'Vapour mole frac.’ (Table X). 

Then Next... 

Observation Tab : Fill in Temperature, Pressure, Liquid density, GOR (Table IX). 

Liquid mol. frac. and Vapor mol. frac. Tabs: Fill in tables with data in Table X. 



Adding experiments... 


h. Add the experiment Two stages separator (SEPS2): Edit / Experiments... /Add /Separators... then fill in 
the items: Follow the same procedure but defining two stages... 

Observation Tab : Temperature, Pressure, Liquid density, GOR (Table XI). 

Component Tab : Liquid mole fraction and Vapour mole fraction for the two stages (Table XII). 

i. Add the experiment Differential liberation (DL1): Edit / Experiments... / Add / Pressure Depletion/ 
Differential Liberation... then fill in the items (procedure is similar to CCE & SEP: 

General Tab : Zl, specify fluid type, temperature. 

Observation Tab: Pressure, Liquid density. Oil relative volume, GOR, Vapour Z-factor (Table XIII). 
Component Tab : Vapour mole fraction (Tables XIV and XV). 

j. Add the viscosity data using the experiment Constant Composition expansion (CCE2): Edit/ 
Experiments... /Add / Pressure Depletion / Constant Composition Expansion... then fill in the items: 

General Tab : Zl, specify fluid type, temperature (Table XVI). 

Observation Tab: Pressure, Liquid viscosity (Table XVI). 

k. Add the experiment Saturation Pressure (PSAT1): Edit / Experiments... /Add /Single Point / Saturation 
Pressure... then fill in the items: 

General Tab: Zl, temperature (Table I). 

Observation Tab: Saturation pressure (Table I). 



Zl is the global composition of reservoir fluid in Table X 
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Adding experiments... 


Once you have added all the experiments, 
your data tree should look like the Figure: 


Save your file!! !! 


To plot the results of the experiments, right click on the experiment 
you want to plot and select Plot. 



Pi* Edt Vftw Run UlfcK“. WhJw Hcfc 


Or- y M&m % » W ft 51 
3 &\& t y U £1 ^ £ 



DLl Oil rtf vtf 

■ ■ CMcKUrt 

DO WHMl 



CCE2 Liqu id vise 



Main windows for fluid characterization... 



Experiments: Go to Edit / Experiment... A window opens where you can either add a new experiment or edit the existing ones. Another 
alternative to Edit Experiments is to right click on the experiment in the data tree and to select Edit.... For example, if you edit un experiment 
you can observe the temperature and various pressures for which the observations were entered... 

Observations: To edit the observations, go to Edit / Observations.... You can practice this point after performing the "Adding experiments" 
instructions! You open a window, which will appear as shown in the figure. On the left the Experiment Type is defined. 

The ones that have already been defined for the current fluid model and that appear in the data tree are marked by an asterisk (*). 

When selecting one of them, the next Experiment list will define the experiments for each category. For example, if you have several 

experiments on the data tree: When selecting one, 


the list of corresponding observations appears. 

The ones marked with * are already selected and entered. 
More can be added or some can be removed by using 
the + and - buttons of the menu of the Observations panel. 
By selecting them, you can see the experimental data that 
have already been entered. You can add an observation 
for viewing even if there is no corresponding experimental 
data available. 

Close the window when you finish looking at the 
observations! 
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Running a simulation 


pi 


Go to Run / Simulate or alternatively, click on the green Go button . You can compare the 
observed and calculated observations for each experiment defined in the data tree. For example, 
compare the observed and calculated saturation pressures at the end. 


Grouping pure components 


Grouping will help to match the laboratory data and obtain a better fluid model. While grouping, 
we need to ensure that the phase envelope is not going to be much affected. 

a. Go to Edit/ Fluid Model / Group.... We will group by mole fraction as selected by default. For 
the components we will group together, we will enter an index in the "New Index" column. 
Enter 1 for IC4 and NC4, 2 for IC5 and NC5, and 3 for C7, C8, C9, and CIO. Select Update. A 
default group name appears in the Group column. Select OK to exit the window. 

b. Check on the fluid composition. You now have 10 different components. When checking in 
the components window (Edit / Fluid Model / Components), you can see that the 
components you have just defined are now of the type User defined, with their 
characteristics in the User tab. 
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Grouping pure components 

c. We will now compare the phase envelope to the one before grouping. 
Select Superimpose and plot the new phase plot in the same graph as the 
previous one. You can see that the two phase plots match very well. You 
can also superimpose the two fingerprint plots to review your grouping. 

d. Save this project as Grouped. PVI (in the menu File /Save os...). We will use 
the fluid model with the grouped pseudo-component for the regression. 


Performing a regression 

a. Use the data-set you have just saved after grouping Grouped. PVI . 

b. We will first set the weight of the regression for all the observations: 

• Saturation pressure is a very important criterion for an EOS match! So 
set the weight of Psat as 60 (right-click on Sat. pressure). 

• Viscosity is always the last to be fitted in a regression... So set its 
weight to zero, right-click on CCE2 and select Don't use in regression. 

• Set the weights of the other observations as follows: CCE1 relative vol. 
to 20, Liquid density to 50; SEPS1 Gas-Oil Ratio to 5, Liquid density to 
5; SEPS2 Gas-Oil Ratio to 5, liquid density to 5; DL1 Oil rel. Vol. to 40, 
Liquid density to 15, Gas-Oil ratio to 15, and vapor z-factor to 15. 
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Performing a regression 


c. 


d. 


Before you proceed with the regression, plot the different experiments and 
select the ones you would like to monitor (there is no space for all of them). 
A good selection is to monitor the CCE1 experiments, the SEPS2 Gas-Oil 
ratio and the DL1 experiments. 

Go to Run / Regression... Select the Special regression (MW of characterized 
components etc...). Select variable and enter 1 for the C11+ MW regression, 
OK, Run! 


H Rtgr«si o n Psrtd 
VwaWw 

Tvp«dii9n»i vanjfWw 


II a S^l 


fomtat (cQTTwncrt pegpetei 6Cs} 




ftp 


Select Special EOS parameters for reg~. |e£jnl 

f^olecUar Wetgf** of Charactersed compcnerts 


Observe the curves during the real time regression, the fit of the experimental data is sometimes 
improved. Check the regression report by clicking on Regression in the Regression Panel window. 
Check the new values of the calculated saturation pressure and all other observations. Notice that 
DL1 oil relative volume is improved but still needs to be worked on. Go to the other tabs of the 
Regression Report Panel. In Modifiers, observe the change that was made on the mole weight of 
C11+ to achieve this fit. It is very important to check the modifiers when doing a regression, to 
ensure no unrealistic or too important changes were made. Accept the regression! 

Save the project as Special. PVI ! 



MW 

+ 

C11 + 

1 

— 
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OK 


Cancel 


’Enter Group Indies For MolectJar WetgNs 




Performing a regression 


g. We are now going to improve the regression on file Special.PVI using other fitting parameters. In the "Regression Panel" 
select the Normal (component properties and BICs) regression variable type. 



h. Critical parameters are well defined for pure components, so we should avoid selecting them. In general, when deciding 
which variables to choose, consider the following sensitivities: critical pressure and temperature of Q a and Q b for the 
saturation pressure; Volume shifts for Z-factor, densities, etc. The rules for the matching parameters selection are as 
follows; In the same column , entering a same number for different parameters means they will be varied dependently. A 
different number means the parameters are matched independently. For the top table of the variable panel , all columns 
are independent , even if the same number is selected in the same row. For the bottom table (BIC variables ), number must 
be different in a row for the parameters to be matched independently. 

i. We usually focus at first on pseudo-components. 

Use the values shown in the Figure... 


Regression Panel Gs H s @ 

vurffcs 

Type of regresMn vanabtes 
^ Normal (component properties and GlCs} 

Special (MW of components etc } 

> PVTt selects 

Variables [ Unto 

Report 

Regression | Smuldbon [ 

Regress 

Rm Accept | Reject | 

Pose | Help | 



Performing a regression 

g. Close the window and click on Run to run the regression! 



h. 


h. 


Observe the calculated versus observed results on the graphs. Check the regression report (Regression button in the 
Regression panel). Most observations are now matched to an acceptable limit. Density is within less than 1%. The 
saturation pressure is not as matched as it was after the Normal regression, but it is still within an acceptable range. 
Accept the regression. 


We now need to match the viscosity data. We first remove all 
the observations from the regression using the Don't use in 
regression by right-clicking on each experiment (CCE1, SEPS1, 
SEPS2, DL1 and PSAT1) except that of the viscosity 
observation. Set the weight of the CCE2 liquid viscosity to 1, 
making sure that CCE2 is Used in regression. Plot the CCE2 
liquid viscosity to follow the regression in real time. Go to the 
Regression Panel and select Variables. Clear all the selected 
regression variables by clicking Reset. Enter values as in the 
figure for VcritV (C6 and C7+ have the same number to make 
them dependent during the regression). 


Run the simulation. Check the match from the Regression 
Report Panel. 
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Viscosity match... 




Performing a regression 



NOTE: 


PVTi 




If you get a warning message informing you that the model is non-monotonic. Check the log area for 
the reason. The message will inform you that the way the variables were varied, C7+ critical volume for 
viscosity is less than that of C6. Reject the regression. 




The flwd model ■$ not moronic 
Seethe Log WWtdowtor rformebori. 

Failure to meet mooof&nictty eonatrainls js 
usualy an ndcauon that the model w® not be a 


umIuI fluid descnptfon. 


In the Limits panel put all the regression limits to between 0.5 and 1.5. Run the simulation. Check the 
regression report. This match is even better than the previous one: plot the experiments!!! 



Accept the regression. 


Save your project as 'Matched. PVr 
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Exporting water correlations 

From your Matched. PVI file: 

a. Go to File/ Export Keywords / Water.... Enter the reservoir temperature of 111 ° C and the reservoir pressure of 445.0 
bar(g). The water salinity is around 17,000 ppm in this area. Because the exact value is not given accurately in the PVT 
report and because this amount of salt will not change the output much, we will not enter further information. 



U PVTi 

Y*t\ E<tt Wp* 



Uliditief Qptwni Window Mdp 

kjz w 

uw<* c PS 


ECllPSE CflmiHJsitioflil fluid Model-. 
ECUKE Thermo Flurf Model... 

APt Tricking option n EtUPSE tl* *QiL. 
Oil rtiefv-ou. .. 

C«s W«vgii_ 

Equilibrition-., 

Wdit., 

YfPi 



b. Click OK and save the file as PVTW.PVO when prompted. The information you have just saved is displayed in the Output 
display window. PVTW is the keyword that will be used later in the ECLIPSE simulation for the water PVT properties. The 
different arguments of the PVTW keyword are the reference pressure, Bw, Cw and pw at reference pressure, and the 
"viscosibility" = (dpw/dP)/pw) 
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Exporting ECLIPSE BO tables 



We will now export the ECLIPSE Blackoil keywords of the fluid model we have just matched, as tables. From your Matched. PVI 
file: 

a. Right click on DL1 and select Export Keywords... Select PVTO & PVDG (Live oil and dry gas) properties keywords. Select the 
separator Experiment SEPS2 and click OK. Save the file as PVTO.PVO when prompted. 





b. The ECLIPSE Blackoil keywords are now created and appear in the Output Display window. Export keywords for ECLIPSE 
Blackoil are computed using the DL1 pressures. 

c. Look in the Output display window. The keywords: 

DENSITY provides the densities of 0, W and G at surface conditions. 

PVTO (live oil with dissolved gas) provides the dissolved gas (Rs), pressure. Bo, and p. As you can observe, 
the table consists of different sub-tables, each one for a given Rs corresponding to the Psat defined in the 
experience DL2. 

PVDG provides the pressure, Bg and pg (the pressure are the values provided in the DL1 experiments). 
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Exporting ECLIPSE BO equilibration keywords 


We now need to export the keyword for the definition of Rs versus depth. To generate the equilibration keyword, we need 
to define an experiment of composition versus depth. 

a. Go to Edit / Experiments and Add Composition with depth.... The fluid Zl is already selected. Click Next. Change the 
fluid type to Oil. Enter the reservoir temperature 111 ° C. Enter the reservoir pressure as 445 bar(g) and the depth as 
3160 m. In the observation tab, enter the depth 3100 m in the first row and 3240 m in the second. The intermediate 
depth is automatically generated and filled by ECLIPSE PVTi. Click Next. Say OK to the message, the experiment 
COMPG1 is now included in the data tree! 
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b. Right-click on this new experiment and select Export Keyword.... Select PBVD/PDVD equilibration keyword. Save the file 
as PBVD.PVO when prompted. 

c. Look in the Output display window. The keywords: 

PBVD (Bubble point versus depth) provides the depth entered in the COMPG1 experiment and the bubble 
point pressure. 
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Dynamic reservoir simulation 
Data review: Initialization 



Data review: initialization 

► Initial state 300 

► Contacts 310 
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Initial state: summary 



► The reservoir is in static equilibrium 

► The objective is to set up: 

• Initial pressures 

• Initial saturations 

► The initial state is correct when original volumes in place are correctly 
calculated 

► Main data influencing Original volumes in place are: 

• Net Pore volume 

• Saturation Height Functions (influenced by Pc curves and fluid densities) 
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Initial state: water-oil contact definition 



► Three main intervals can be defined: 

• HC leg (Sw = Swc = Irreducible water saturation) 

• Transition zone (Swc < Sw < 1) 



Capillary pressure curves 



The model initial equilibrium saturation 
distribution is established with drainage 
capillary pressure curves. 

Swi: irreducible water saturation 

(* inital water saturation) 

Sor: residual oil saturation 
Pcd: Displacement cap. Pressure 



0.0 0.2 0.4 0.6 0.8 1.0 


Water saturation, S w _ 1 ■> Oil 

w Sor 
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Capillary pressure in the reservoir 



Oil only 


Oil + water 
Water-only 



Pc 


Well 


Pc = Rho x g x h 


RFT 

gradients 


Pc^ h 


Swi = f (h) 


Sor = f (pc) 


Transition 

zone 


WOC 


Capillary 

pressure 


Reservoir 


FWL = Free Water Level -> Pc = 0 

WOC = Water-Oil Contact Displacement pressure (Pc d) 


RFT / MDT FWL 
LLD / MSFL -> WOC 
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Initial state: pressure calculations 



► Pressures gradients are proportional to fluid densities 

GOC 

WOC 

► Pressures changes at interface can be derived form Pc curves 

• At the gas oil interface: Po (GOC) = Pg (GOC) - Pcgo (GOC) 

• At the water oil interface: Po (WOC) = Pw (WOC) + Pcwo (WOC) 

► In order to calculate the initial pressure, it is sufficient 

• To give the pressure at a reference depth 

• The location of the contacts and the capillary pressures at the contacts 
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Initial saturation distribution 



Po = Po ref + p 0 gZ 


Pw = Pw ref + p w gZ 


Po — Pw 


Sw in the transition zone 
OOIP calculation 
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Initial pressure distribution 



Hypothesis: 

► Reservoir at initial equilibrium (most often gravity): 


ll 



Po Po re f + p 0 gZ 


Po known everywhere 


► At WFL: Pc = 0 and then Pw FWL = Po FWL 


Pw known everywhere 
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Pw = Pw ref - p w gZ 






Initial state: saturation calculations 




► There is a gravity equilibrium in the transition zone 

• In the water leg: S w = 1 

• In the oil leg: S 0 = S 0 (Pc wo ) with Pc wo = P 0 (z) - P w (z) 

• In the gas leg: S g = S g (Pc g0 ) with Pc g0 = P g (z) - P 0 (z) 

► Saturations are calculated at the grid block center and are applied to the 
entire grid block 
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Example of the initial reservoir condition calculations 


► Goal: 

• Compute pressure and saturation distribution everywhere within the reservoir model at 
the simulation start 

► Variables: pressure (P) and saturations (S) 

► Example: 

• Simple two-phase flow case oil/water 

• Initial data: 

• WOC 

• P @ reference depth z (datum) 

• Compute pressure distribution 

• Compute saturations 
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Initial state 



HAFWL 
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Contacts discretization 





Exercise 


Given the following data for a reservoir (single rock type) containing oil and water 
only:| 

1 Calculate the height of the water-oil transition zone, assuming constant reservoir 
oil composition. 



Pressure at datum, Pi = 59 bar 
Bubble point pressure, Pb = 57 bar 
Bo (Pi) = 1.1049 Bo (Pb) =1.105 
Surface Densities: p 0 = 896 kg/m 3 


Bw (Pi) =1.03 Rs = 32 v/v 
p w = 1059 kg/m 3 p g = 0.9 kg/m 3 


Sw 

Pcow (bar) 
(drainage) 

Pcow (bar) 
(imbibition) 

0,11 

0,6 

0,6 

0,20 

0,2 

0,15 

0,30 

0,1 

0,04 

0,40 

0,09 

0,0 

0,50 

0,08 

-0,04 

0,60 

0,07 

-0,1 

0,70 

0,06 

-0,5 

0,77 

0,03 

-1,0 

1,0 

0 



2. Calculate the height above the free water level at which the oil saturation is 50%. 
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Initial distribution 


► Several traps (non-communicating zones of the reservoir) can be defined. The 
use of several traps allows to take into account: 

• distinct GOC and/or WOC 

• distinct Pb vs depth laws 

• distinct initial pressures 

► In the case of PVT traps: 

• A complete PVT model (oil, gas and water characteristics) is affected 

• Pb vs depth curve, GOC, WOC and pressure at a given depth (datum) are given 

NB1 : initial pressure curve, Pb curve and GOC must be compatible 

NB2 : knowing GOC and WOC depths and fluid densities, the simulator can 
calculate initial pressure vs depth curve. 
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Initial state with several traps 
Multi layer reservoir example 



► A trap is defined as a set of cells in static equilibrium 

► A trap is described by: 

• Water-oil contact and a value of Pc wo (WOC) 

• A gas-oil contact and a value of Pc g0 (GOC) 

• An initial pressure at datum 

► A flow model can be initialized with one or several traps 



Key points to keep in mind 



\ 

► The initial state is correct when original volumes in place are correctly calculated 

► Three main intervals can be defined: 

• HC leg (Sw = Swc = Irreducible water saturation) 

• Transition zone (Swc < Sw < 1) 

• Water leg (Sw = 1) 

► The model initial equilibrium saturation distribution is established with drainage 
capillary pressure curves. 

► Pressures changes at interface can be derived form Pc curves 

• At the gas oil interface: Po (GOC) = Pg (GOC) - Pcgo (GOC) 

• At the water oil interface: Po (WOC) = Pw (WOC) + Pcwo (WOC) 

► In order to calculate the initial pressure, it is sufficient 

• To give the pressure at a reference depth 

• The location of the contacts and the capillary pressures at the contacts 

► Several traps (non-communicating zones of the reservoir) can be defined. 
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Dynamic reservoir simulation 
Data review: Aquifers 


IFPT raining 



Data review: aquifers 


► Aquifer description 320 

► Aquifer models 330 
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Aquifer 

description 


West edge 
Aquifer 


Aquifer representation 


Objectives 

► to know cumulative water influx and pressure support versus time 

► to approximate water influx by a model 



Aquifer representation 


► Objective is to estimate the aquifer activity, i.e. to estimate cumulative water 
influx and pressure support versus time 

► These results depend on the following parameters: 

• Vaq = Aquifer Volume 

• Ct = Aquifer Compressibility 

• Taq = Transmissivity between reservoir and aquifer 

► There are 3 possibilities to describe an aquifer: 

1. Describe the entire aquifer with a large number of cells 

2. Use large grid cells 

3 . Use Hurst & van Everdingen functions 
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Aquifer's configurations 


► The bottom water drive aquifer: it is in contact with the entire hydrocarbon 
area. Water invasion occurring vertically is governed by the reservoir vertical 
permeability. 



► The edge water drive aquifer: the water entries take place laterally. 
Horizontal permeability is governing the water mouvement. 



WELL 


AQUIFER 


► The linear water drive aquifer: it is a typical shape encountered 
in fluvio-deltaic deposition reservoirs (channels) 



Oil 


Aquifer's representation 



► Water influx arises as a result of pressure changes decompressing the aquifer. 
If the pressure in an aquifer can be calculated, the resulting volumetric 
changes can be determined from the pressure/volume compressibility 
relationship. 

► Pressure profile in a reservoir - aquifer system: 


Pressure 
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Aquifer's representation 



► Cumulative water influx and pressure support versus time depends on the 
following parameters: 

• Aquifer volume: Vaq 

• Aquifer compressibility: Caq 

• Transmissivity between reservoir and aquifer: Taq 

► There are two possibilities to represent aquifers 

• Gridding aquifers 

• Numerical aquifers 

• Analytical aquifers 
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Main parameters 



► Aquifer expansion: 


► Aquifer compressibility: 

• Cp: pore compressibility 

• Cw: water compressibility 



► Aquifer depletion: 
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Aquifer's modeling 
The diffusivity equation 



► The combination of the mass conservation law with Darcy law lead to the 
diffusivity equation: 

f d 2 P ldP^i , ^ udP 
^ dr 2 r dr J k dt 

> Assumptions: 

• Radial flow 

• Homogeneous isotropic porous medium, uniform thickness of the medium 

• Constant porosity and permeability (independent of pressure) 

• Fluid of small and constant compressibility, constant fluid viscosity 

• Small pressure gradient, gravity neglected (radial axis) 

THE SOLUTION OF THE DIFFUSIVITY EQUATION DEPENDS ON INITIAL AND BOUNDARY CONDITIONS 
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Aquifer's modeling 

Initial and boundary conditions 


► Initial condition: P = Pi @ t = 0 everywhere in the aquifer 

► Inner boundary condition: the most commonly used, 

• the constant terminal rate, P = f(t) i.e. used to establish build-up analysis 

• the constant terminal pressure, q = f(t) i.e. used to calculate water influx 


► Outer boundary condition: 

• infinite aquifer, P = Pi as r -> 00 


bounded aquifer, 





Aquifer models 




Aquifer Representation 




The usual boundary condition to grid is "No flow" 


► Three possibilities to describethe aquifer: 

1. Describe the entire aquifer with a large number of cells 

2. Use large grid cells 

3. Use Hurst & van Everdingen functions. 



Aquifer representation 



► 

► 


The usual boundary condition to grid is No-flow 

To represent the aquifer, the following possibilities exist: 

• Grid the whole system (possibly with larger cells in the water zone, where lower 
precision is needed) 

• Take into account only a limited number of water cells in the WOC vicinity (one or a few 
rows) and multiply their pore volume (and correct their permeability if necessary) 

• Take into account only a limited number of water cells in the WOC vicinity (one or a few 
rows) and implement an analytical or numerical aquifer in these cells 

The last two options allow for reducing drastically the total number of cells 

(they are usually the two options between which the engineer has to make a 

choice) 
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Use of large grid cells 



N 


Or a large number of grid cells 



1 km 

I 1 
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Numerical aquifers 


► Instantaneous expansion of a finite volume 

► This is usually applied to small aquifers with closed boundaries and in 
equilibrium with the reservoirs. 

• Qwi = CdP/dt 

• C = dV/dP 

• This is equivalent to a multiplied pore volume 
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Analytical aquifers 
Pseudo-steady state: Fetkovitch 



► The Fetkovich aquifer model uses a simplified approach based on a pseudo- 
steady state productivity index and a material balance relationship between 
the aquifer pressure and the cumulative influx. 

► This model assumes that the pressure response is felt uniformly throughout 
the entire aquifer: the pressure at the outer boundary of the aquifer does not 
change. 

► The aquifer inflow is modeled by the equation: 


Q wi = C(Pi-P) 

C: equivalent to a transmissibility 
Pi: initial pressure in the gridblock i 
P: pressure at time t in the gridblock i 




= P 


INI 


P 


RES 


* 


r 


IN 


► 
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Analytical aquifers - Transient Aquifer 



► All aquifers are finite in size, however there is a period of time when a 
pressure disturbance created by production from a well has not traveled far 
enough and reached the boundary of the aquifer. 

► During this time the aquifer behaves as being infinite and unsteady state flow 
applies. After the boundary influences the behavior of the system and 
pseudo-steady state flow starts. 
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Analytical aquifers 
Transient aquifer: parameters 

► aquifer geometry 

• netthickness 

• aquifer inner radius (radius equivalent to contact length) 

• ratio r EXTERNAL / r, NNER = ro/ri 

• encroachment angle (aperture) 

► petrophysics 

• porosity 

• permeability 

• compressibility of water bearing rock 

► water 

• compressibility of aquifer water 

• viscosity of aquifer water 

► external boundary condition 

• condition at ro (constant pressure or no flow) 
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Radial aquifer model 
Hurst & Van Everdingen 

Hurst and Van Everdingen solution is 
applicable to any reservoir where the flow 
into the reservoir is essentially radial in 
nature. 

Hurst & Van Everdingen solutions are 
functions of: 

• Dimensionless water influx and 
dimensionless time, as a function of 
dimensionless radius. 

In well testing, the diffusivity equation is 
solved between the wellbore radius and the 
drainage radius, the same method will be 
used, but between the WOC radius and the 
external aquifer radius 

!FP Training I 338 


Analytical aquifers 

Transient aquifer: Carter-Tracy model 


► The Carter-Tracy aquifer model is a simplified approximation to a fully 
transient model, which refers to the diffusivity equation. 


► The method uses tables that supplie a constant rate and a constant pressure 
influence function: Hurst & van Everdingen. 

► Although the theory has been developed for a radially symmetric reservoir 
surrounded by an annular aquifer, the method is applicable to arbitrarily- 
shaped reservoirs. 
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Analytical aquifers 

Carter-Tracy model: Boundary conditions 



Two cases: 

► The Pressure case, where the pressure at the inside boundary is known and 
there is no flow at the outside boundary and we calculate the water influx. 


► The Rate case, where the rate is known at the inside boundary. At the outside 
boundary there is no flow and we calculate the total pressure drop. 
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Analytical aquifers 

Carter Tracy model: parameters 





Water influx: We = Urjr [y 6, k.h, </>.h.C) 

r 0 : outer radius k\ permeability 

r,: inner radius </r. porosity 

O encroachment angle h : thickness 

C: compressibility 
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Analytical aquifers 



Carte Tracy: dimensionless parameters 

► The two main parameters that govern the behavior of the aquifer are: 


the dimensionless time: 



k 

f 

|iW 

C 

ri 


the aquifer influx constant: 


U = h-0-(/)-C-r: 


h 

q 


aquifer permeability 

aquifer porosity 

water viscosity in the aquifer 

compressibility 

outer radius of the reservoIR 


aquifer thickness 
aperture 
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Analytical aquifers 



Carte Tracy: constant pressure case 


P*<R.O) = 1 


0<R.t*> 

~ a R = 00 R=3.0 



*)=0 


Q R (t*)=We(t)/2x 

outer radius: R permeability: 1 

inner radius: 1 porosity: 1 

aperture: 2n thickness: 1 
compressibility: 1 



1 2 3 4 5 6 7 


We = e ■ <t> ■ h ■ C ■ r ; 2 • AP • W D ( t * ) 
We = U ■ AP • W D (t*) 
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Analytical aquifers 



Carte Tracy: constant rate case 



*')=!- AP*{t*) 



outer radius: R permeability: 1 

inner radius: 1 porosity: 1 

aperture: 2 n thickness: 1 
compressibility: 1 


p( D t*} - t*} 4 - ^ - Fi 

1 1 

L )-r 

2-6 -k -h 

4 -t*J 
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Key points to keep in mind 



<Xu 


► Gridding the aquifer? 

• Time and memory consuming 

► Porosity "increase" of the last cell row (or consider a CLOSED aquifer)? 

• All the water is located at the same distance 

• What increase of Pore Volume? 

• How to correct permeability? 

► Introduce STEADY-STATE aquifer? 

• Only applicable to very particular cases corresponding to constant Pressure (close 
outcropping, drastic change in permeability, partially sealing fault...) 

► Use "transient" aquifers? 

• Reproduce the propagation phenomenon 

• Physical parameters (identical to the reservoir ones except for external radius) 
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Dynamic reservoir simulation 
Data review: Flow description 



Sommaire 


► Introduction 348 

► Flow description - one phase 352 

► Flow description - multiphase 370 

► Hysteresis 390 

► Options in reservoir simulation 402 
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Introduction 


Reservoir zones: 

► Mono phase: 

• Oil leg 

• Water leg - aquifer 

► Multiphase 

• Transition zone 


RESERVOIR 


PRODUCTION BEHAVIOR 


CAPILLARY PRESSURE, Pc 



Pc y 


Sw = Sw c 

Clean Oil 

No Water Cut 


1 > Sw > Sw c 

^ Increasing Water Cut 


Sw = 1 

Water Only 



RESERVOIR CONDITIONS 


non -volatile oil 


gas condensate 
huons p 

critical point j 

volatile oil ^ ^ 

a / i 


wet gas 

P 

f 

/ 


Displacement in production 


monophase or multiphase flow 


Reservoir fluid - PVT: 

► Mono phase: 

• Undersaturated oil 

• Wet gas (above dew point) 

• Dry gas 

► Multiphse: 

• Saturated oil 

• Gas condensate (below dew point) 


SEPARATOR CONDITIONS 


TEMPERATURE 




Permeability definition 


Dimensions 

► Measures the capacity and ability of the formation to transmit fluids 




1 cp - _L x dyne_y_sec 


Darcy's Law 

p: Fluid Viscosity (cP) 

AP: Differential Pressure (atm) 

A: Cross sectional area (cm 2 ) 

Q: Injection flow rate (cm 3 /s) 

L: Length (cm) 


cm 3 /sec 
1 Darcy \ 


100 cm 2 


cm 


cm 


A />-/> 

/ \ _ 

1 atm =1.033 kg/cm 2 
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Relative permeability 



► Darcy Law in monophasic fluid flow 


Q = -K 


A 


H • AL 


AP 



► Darcy Law in multiphase fluid flow 


AL 


Qo = 


k.A AP. 


Po AL 


Q* = 


k.„ A AP... 


kr. =— 


AL 


K 


• K = absolute permeability (m 2 for SI, mD for PFU unit systems) 

• ko = effective permeability to oil (same units as K) 

• kro = relative permeability to oil (no dimension) 
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Flow description - 
one phase 




Monophasic permeability 



► Single phase flows are modelled with Darcy 's law: 


U = grad O = 

ja ja 


dP 


V 


dx 


Pg 


dz 

dx 


A 


J 


► Permeability K is independent from the fluid 


► Laboratory measurements are done on plugs which are extracted horizontally 
or vertically from cores. Experience shows that Kv is not equal to Kh 

► Gas and oil flow with connate water in the reservoir 
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Permeability dimension 



[K] = [U] . [*i] / [dP/dx] 

[K] = [L T 1 ] . [M L 1 T 1 ] / [M L ' 2 T' 2 ] 

[K] = [L 2 ] 

Permeability has the dimension of a surface 
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Flow calculations between cells 
Transmissivity 



[ T ab ] = [Q ab ] . [p / [P ab ] 


\ab 


ab 


ij,k+l 


i, j-1, k 


a 




b 


i+l, j, k 


i. j> k 


i,j,k-l 


ij+l,k 


[ T «b] = ^ T -ij [ ML -t r'] /[ML' T' 2 ] 
[T ab ] = /i 3 ; 

Transmissivity has the dimension of a volume. 


Qab = - (T“> / H ) ( AP - p g Az) 


T=K ± 

L 


K : reservoir permeability 
A : surface open to flow 
L : length of flow path 
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From core to flow model 



core 


Grid block 
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Single phase flow: problems to solve 


► Permeability composition 

• Problem: how to represent a permeability distribution by a constant permeability when 
the flow path is known 

• Solution: arithmetic average, harmonic average, power law... 

► Permeability homogenization 

• Problem: how to represent a permeability distribution by a constant permeability when 
the flow path is unknown 

• Solution: fine scale model, history match... 

► Permeability modelling 

• Identification of basic laws (e.g.: Phi - K laws) 

• Well blocking & Upscaling of basic laws 

• Population of cells 

• Quality checks 
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Permeability composition: random isotropic flow 




Q 

Ln K = 2 pi Ln Ki 


Permeability composition: horizontal & vertical flows 
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Vertical flow 
ID flow 

H / Kv = D hj / KVj 


Horizontal flow 
Parallel flow 
Kh H = 2 Kh: h: 
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Permeability composition: power law 



► Average permeability is calculated as follows: 

K"=2> K i" 

► This is an extension of harmonic, geometric & arithmetic means: 

• w = -1: harmonic mean 

• w = 0: geometric mean 

• w = 1: arithmetic mean 

► Coefficient w is related to the type of flow 

• w = -1: serial flow 

• w = 0: random isotropic flow 

• w = 1: parallel flow 
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Permeability composition: power law 



G) = -l 

harmonic mean 
serial flow 


co = 0 

geometric mean 
random flow 


co = 1 

arithmetic mean 
parallel flow 



This is one example of composition, the shape of the curve depends on the 
K distribution 
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Permeability composition: transmissivity calculations 




ab _ 2 Kx a Kx b AyAz 
~Kx a + Kx b Ax 

ab _ 2 Ky a Ky b AxAz 
_ Ky a + Ky b Ay 

ab _ 2 Kz a Kz b AxAy 
_ Kz a + Kz b Az 
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Permeability composition: transmissivity calculations 



Tz 



Tx 


For a regular grid and a homogeneous reservoir: 
Tz / Tx = (Kz / Kx) / (Dz / Dx ) 2 
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Permeability homogenization: principles 




^3 



Q4 


HETEROGENEOUS CELL 


HOMOGENEOUS CELL 
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Permeability homogenization: some results 



► Cardwell & Parsons (1945) 

• Hypothesis 

• K is a random function 

• Results 

• Upper & lower limit for K* derived from K distribution 

► Matheron (1966) 

• Hypothesis 

• K is an isotropic random function 

• Results 

• K* derived from flow dimension, K harmonic & arithmetic means 

► Begg & King 

• Hypothesis 

• Reservoir is composed of sands (constant K) and impermeable shales 

• Results 

• K* derived from K and shale distribution 
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Permeability homogenization: Matheron formula 


► Hypothesis 

• Reservoir is infinite & isotropic 

• Flow dimension (ID, 2D, 3D) is known 

• K distribution is known 

► Results 

• Equivalent permeability depends on flow dimension, K harmonic & arithmetic K means 

N - 1 1 

k* A(K) + — H(K) 

N N 


A(K) = X P i k , 


H(K) = 1/ X 


P_i. 

k ; 


Matheron formula respects the range H(k) < K* < A(K) 
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Key points to keep in mind 


<Xu 



► Single phase flows are modelled with Darcy 's law 

Problems: 

► Permeability composition: how to represent a permeability distribution by 
a constant permeability when the flow path is known 

^ Permeability homogenization: how to represent a permeability 

distribution by a constant permeability when the flow path is unknown 

► Permeability modelling: Identification of basic laws, upscaling of basic 
laws... 
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Flow description - 
multiphase 


Displacement in Production 


► The displacement will depend on: 

• Relative permeabilities 

• Oil and water viscosity 


RESERVOIR 


Sw = Sw 


1 > Sw > Sw 


PRODUCTION BEHAVIOR 


Clean Oil 
No Water Cut 


Increasing Water Cut 


Water Only 


CAPILLARY PRESSURE, Pc 


Oil and water flow at equal ease. 

Residual oil saturation (Sor), this oil cannot 
be recovered by water flooding. 

Non movable water saturation (Swc) 
characteristic for a specific rock. 


o Sw, 


Sw=l 


Sw = 1-So 
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Multiphase flow description 



► Polyphasic flows are modeled with the following equations: 

• Ug = -K(Krg/|ig)grad®g 

• Uo = - K (Kro / \xo) grad <Do 

• Uw = - K (Krw / |iw) grad ®w 

► Relative permeability curves measured in the laboratory are based on ID - 
two phase flow (water/oil or gas/oil) 
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Multiphase flow 
Relative permeability curves 



Kr Kr 




WATER-OIL GAS-OIL 
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Multiphase flow 

Relative permeability normalization 



Kr KrD 




Sw - Swc Krw 

SwD = KrwD = 

1-Swc-Sorw Krwmax 


Kro 

KroD = 

Kromax 
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Multiphase flow 

Relative permeability normalization 



KrD 



2 SwD 


Kr 



Sw = SwD (1-Swc-Sorw) + Swc 
Kro = KroD . Kromax 
Krw = KrwD. Krwmax 
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Multiphase flow 
Calculations between cells 



I, J, K+l 



► Monophasic flow: Q ab = - (T ab / p) . (AP - p g Az) 

► Polyphasic flow: Qj ab = - (T ab Krj ab / p ; ) . (AP S - g Az) 
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Multiphase flow calculations 
Use of upstream saturation 



1. Pressures and saturations are defined at the grid cell centre 

2. Flows are calculated per phase: 

• Q q = -(k. A /Ax) . M g . (AP g -p g .g. Az) 

• 0 o = -(k. A /Ax) . M 0 . (AP 0 -p 0 .g. Az) 

• Q w = -(k. A /Ax) . M w . (AP w -p w .g. Az) 

3. For each phase, the flow direction is defined by the potential gradient 
between two grid cell centres. The phase mobility is calculated using the 
upstream saturation 
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Multiphase flow 
Problems to solve 

► Three phase relative permabilities 

• How to represent a three phase flow with two phase rel perms 

► Homogenization 

• How to represent an heterogeneous porous medium by an equivalent homogeneous 
porous medium. 

► Pseudoization 

• How to represent a fine gridded homogeneous porous medium by an equivalent coarse 
gridded homogeneous porous medium. 
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Multiphase flow 

Three phase relative permeabilities 



► Three phase flow are not usually performed in laboratory conditions 


► Three phase relative permeability are calculated using water-oil relative 
permeability (Krw,Krow) and gas-oil relative permeability (Krg,Krog) 


► It is generally admitted: 

• Three phase Krw = Krw (Sw) 

• Three phase Krg = Krg (Sg) 

• Three phase Kro = Kro ( Sw,Sg ) 
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Three phase flow: ternary diagram 





!PP Training 


380 


Three phase relative permeabilities: channel flow theory 



► There is at most only one mobile fluid in any flow channel 

► The wetting phase is located primarily in the small pore spaces, the non 
wetting phase in the large pore spaces, and the intermediate phase spatially 
separates them 

► At equal water saturations, the microscopic fluid distributions at the water oil 
interface are identical in a water-oil system and in a three phase system: Krw 
and Pcwo are identical in the three phase and two phase systems 

► At equal gas saturations, the microscopic fluid distributions at the gas oil 
interface are identical in a gas-oil system and in a three phase system: Krg 
and Pcgo are identical in three phase and two phase systems 
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Three phase relative permeability: Stone 1 model 



kr a = 


S c - kr oW kr o g 


kr 


owmax 


a-sj.a-s g ) 


► Normalized fluid saturations: 


£* _ S W C 


(i - S wc - S om ) 


£* _ S 0 Som 


(1 - - S om ) 


O* _g 

g (l - s - s ) 

\ wc om / 


S om = minimum value of the residual oil saturation 

Q _Q 

Q -Q orw ^org o 


y om ^orw 


i_« _s g 

^wi org 


► Impedance factors of oil flow: 

Kr„=fi ow {l~S:) Kr og =P og {\-S' g ) 

► Impedance of oil flow by water and gas are two independent events: 

Kr„ = 
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Three phase relative permeability 


Ternary diagram for isoperm representation (Stone I) 



Sg * 100% 


Som 


Sw= 100% 



Vs^-ioo% 
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Three phase relative permeability: Stone 1 (example) 





Geometrical construction 


o 



W G 


Swc = 20% Sgc = 5% 

Sorw = 20% Sorg = 15% 
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Three phase relative permeability: Stone II 



(kr w + kr g ) 

► krow(Sw) and krog(Sg) are computed from the 2 phase flow kr since the third 
fluid is absent So = 1-Sw or So = 1-Swi-Sg 

► Som is not computed 

► kro may be negative, in this case kro = 0 


kr G = max(0, kr ( 


f 


ow max 


kr. 


\ 


ow 


w 

V ow max 


+ kr, 


w 


f 


J 


kr. 


A 


og 


+ kr n 


IcT* 

ogmax j 
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Three phase relative permeability: Other models 



► Stone-Dietrich model: 


kr o =~jj~ — [ kr oS s J+K^ h '„ s (. s g )+ kr g Y\ kr ™ + k\ 


ow max 

identical to Stone II if krowmax = 1 


► Baker model I: 

weighted interpolation with the saturations 

kr S g k r JS c ) + (S w -SJk ro JS 0 ) 

S s +S w -S wi 

krow and krog are computed with So=l-Sw-Sg 
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Three phase relative permeability 
Comparison 



► Stone's model II underestimates kro and leads to non physical solutions for 
small values of So 

• trapping of oil: Som > max(Sorw, Sorg) 

• or, on the contrary, Som = 0 (mainly if pseudo) 


► Stone's model I avoids the non physical solutions but needs the knowledge of 
Som. Overestimation of kro for small values of So 


► Baker's model I generates non linear isoperm if Sorw and Sorg are very 
different 


► All models except the geometrical construction diverge for the small So. The 
final recovery is biased. 
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► Relative permeability 

• Allows to generalize Darcy's law to multiphasic flows in the reservoir 

• Is defined by a curve as a function of saturation, typically k ro (S w ) and k rw (S w ) 


► Relative permeability curves 

• Depend on wettability 

• Rule flow of fluids (water and hydrocarbons) in the reservoir and have a strong 
influence on the Recovery Factor 

• Are measure in laboratory in ID and two-phase flow 

► Three phase relative permeabilities are calculated through models: 

• Stone 1, Stone 2 

• Stone-Dietrich 

• Baker 

• Etc... 

► Models are based on normalized fluid saturations, water-oil relative 
permeability (Krw,Krow) and gas-oil relative permeability (Krg,Krog) 

v 
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Hysteresis 




The hysteresis effect 



► Model used to simulate reservoir fluid flow when a non monotonous 
saturation variation exists during production (drainage, imbibition cycles) 

• Important W/O transition zone or strong capillary effect in dynamic 

• Gas storage in aquifer reservoir 

• Gas injection after water sweeping, tertiary recovery process (WAG) 

• Variation of the WOC or GOC 

• Water injection in depleted reservoir having mobile free gas 
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Definition of the hysteresis effect 



► In both water-wet and mixed-wet conditions, hysteresis in relative 

permeability and capillary pressure accompanies changes in saturation. This 
reflects the difference between water-advancing and water-receding contact 
angles, and the location of oil and water in the pore spaces. 
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Definition of the hysteresis effect 




Wael Abdallah et al., 

Oilfield Review, SLB. Summer 2007 


► The primary drainage (red) and 
imbibition (black) curves bound 
the capillary pressure behavior. 
If the direction of saturation 
change is reversed at an 
intermediate saturation. Pc will 
follow an intermediate path 
(green). Another reversal will 
take it back to the drainage 
curve (yellow). 

► This behavior could occur in the 
middle of a transition zone, or 
as a result of oil banking during 
a water flood. 
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Hysteresis modeling 




► Hysteresis between the primary drainage (red) and imbibition (black) curves 
can be represented by a series of scanning curves (gold). Each scanning curve 
represents a different starting saturation point on the drainage or imbibition 
curve, which would correspond to different heights in the transition zone. 
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Hysteresis modeling 



► The user provides the simulator with imbibition and drainage curves 

► The simulator computes the scanning curves depending of the saturation 
history 

• SN: non wetting fluid saturation (e.g.: So) 

• SNrmax: maximal residual non wetting fluid saturation (e.g.: Sorw) 

• SNmax : maximal non wetting saturation (e.g.: 1-Swi with WO) 

• SNrhyst: non wetting saturation when the direction of saturation change is reversed 
(e.g.: drainage to imbibition) 

► The scanning curves are always reversed 

► Several available hysteresis models: 

• Killough 

• Carlson 

• Jargon 
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Pc hysteresis modeling 




■ Should another reversal in the direction of 
saturation change occur while Pc is in a 
scanning state, a new scanning curve is 
generated that scans back to the point where 
the bounding curve was left originally. 


■ New interpolating function is used that 
depends on the new reversal point, the 
remembered reversal point, the present 
saturation and the regression parameter 


► SNr is computed by Land formula: 



nhyst 

^ N 


1 + C.S 


hyst 

N 


C = 


1 1 


rrmax 

^ Nr 


S 


max 

N 
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Pc curves: Recommendations for use in a simulator 



► Use the drainage Pc in static and dynamic 

• Only if the transition zone is small so that the effect is negligible in dynamic (e.g. Pc 
imbibition is « 0 or small layer thickness or weakly heterogeneous model) 

• Interest 

• the data is simpler 

• CPU time is lower 

► Use the hysteresis Pc in other cases 

• input the imbibition Pc curves 

• if no measurement, estimate imbibition curves going from Pcmax to 0 
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kr hysteresis modeling (e.g. Killough) 



► Wetting phase 



□ Wetting phase is treated differently 

□ Computed Imbibition Relative 
Permeability curves are reversible 

□ Killough model predicts the maximum 
wetting relative permeability by 
interpolating between the "known" 
imbibition relative permeability at the 
maximum residual non wetting 
saturation and the value of the drainage 
wetting relative permeability at the 
maximum wetting-phase saturation 
(trapped non wetting saturation equal to 
zero or Snc) 
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kr hysteresis modeling (e.g. Killough) 



► Non-wetting phase 



□ Parametric interpolation for the scanning imbibition curves 
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Hysteresis modeling: remarks 



► Carlson Pc model may result in errors when computing the recovery by 
spontaneous imbibition in fractured reservoir. 

► In Eclipse, capillary pressure hysteresis is always applied using the Killough 
method 

► Carlson kr model: no hysteresis on the wetting phase, only on the non 
wetting phase. 

► Carlson's method produces a scanning curve that is parallel to the imbibition 
curve. When this method is chosen, it is important to ensure that the 
imbibition curve is always steeper than the drainage curve at the same kr 
value. If this is not the case, the scanning curve could cross to the right of the 
drainage curve, which may produce an unphysical situation (SNr< 0). 
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Options in 

reservoir 

simulation 




Relative permeability 
Endscale option 



Objectives 

► Modify relative permeabilty tables in an easy way, kr tables are normalized 
and remain always the same, only the end-points are changed and kr curves 
are then recalculated. 

► It's a useful option in History Match simulations. 
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Relative permeability 
keywords 



KRO: Maximum oil relative permeability 

KRORW: Oil relative permeability at critical 
water saturation Swcr 

KRWR: Water relative permeability at residual 
oil saturation (1-Sowcr) 

kr 

KRW: Maximum water relative permeability 
SWL: Connate water saturation 
SWCR: Critical water saturation 
SOWCR: Residual oil saturation 
SWU: Maximum water saturation 


KRO KRORW KRWR KRW 
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Relative permeability 
Example 1 


SWL= 0.24 
SWCR= 0.35 




Sw 
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Relative permeability 
Example 2 


SWL= 0.20 
SWCR= 0.25 
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Relative permeability 

► Relative permeabilities 

• When two (or more) fluids flow simultaneously in a rock / in a reservoir, this results in permeability 
reduction for each fluid permeabilities 

• Depends on saturation 

► Relative permeability 

• The ratio of effective permeability to absolute permeability 

• Absolute permeability is the permeability of the rock when totally saturated with one fluid 

• Effective permeability is the permeability of the rock when partially saturated with one fluid 

► Relative permeability 

• Allows to generalize Darcy's law to multiphasic flows in the reservoir 

• Is defined by a curve as a function of saturation, typically k ro (S w ) and k rw (S w ) 

► Relative permeability curves 

• Depend on wettability 

• Rule flow of fluids (water and hydrocarbons) in the reservoir and have a strong influence on the 
Recovery Factor 
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Dynamic reservoir simulation 
Data review: Well representation 



Sommaire 


► Introduction 

410 

► Wells representation 

416 

► Numerical PI 

424 

► Well performance 

434 

► Well constraints 

446 

► Exercise 

456 


Wells' controls 
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Bottom Hole Pressure (BHP) 


BHP 
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Different complex wells 






Cluster well Multidrain or multilateral well 



Re-entry laterals from a 
vertical well 



3D Well 
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Vertical vs. horizontal well 



► For a simulator the only differences are: 




Perforation direction 
Pressure qno|Ho®iwji>0tal:ion 

Well 


Vertical 

Well 
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Wells representation 


► The flow from the formation to the bottom of the well is governed by what is 
known as the INFLOW - OUTFLOW PERFOMANCE relationship of the well. 

► Wells Representation - Geometry 

• Wells radii are very small compared to well cells (~1/1000) 

• Wells are considered as source points in the well cells 

► Inflow Performance 

• Geometry of perforated intervals 

• Inflow performance relationship per perforated interval 

► Outflow Performance 

• Pressure drop through tubing and surface installations 

► Well control 

• Production schedule (oil rate, gas rate, reservoir void age...) 

• Minimum Bottom Hole flowing pressure 

• Minimum Well Head flowing pressure 
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Wells 

representation 




Well's representation 
The problem in reservoir simulation 

► Problems: 

• The well is too small to be represented explicitly, typically rw = 0.001 Dx 

• The well borehole cannot be represented explicitly within the reservoir grid 

• Size and geometry of blocks are not well suited to model physics of flow around the 
wells: 

• radial fluxes 

• large gradients of pressure and saturation 

• The high pressure gradient near the well is not usually resolved 

• a special model is required 
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Well's representation 

The problem in reservoir simulation 



► "Solution": 

• Injection and production wells are represented as source or sink terms 
• Continuity equation: 

V(p.V) + ® = ^ 

dt 


• Well equation: 

® = Qwell = f(PW,Pb) 


• where: 

• Pw is the well bore pressure 

• Pb is the grid block pressure 
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Well's representation 



Wellbore vs perforated grid block 

► A well is represented in block 'i 1 through the production/injection at qp 


( Am\ 


V 


At 


y block i 


£[(' n 'i')-(out ’i') — q p j] 

phase p 


► qp must be related to: 

• block pressure and saturation: Pb and Sb 

• bottom hole pressure (flowing): Pw 

► The pressure - rate relationship determines the boundary conditions: 

• History matching: Q imposed -> Pw 

• Forecasting: Pw imposed -> Q 
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Well's representation 

Well bore vs perforated grid block 

Important: 

► Where the well is perforated 

► Grid block intersected by the well 
trajectory 

► Connections 



Calcul of simulator: 

► Pressure for each grid block, Pb 

► Well bore pressure for each 
connection, Pw 

► Well's BHP at a reference depth, 
Pbh 




Well's representation 

Well bore vs perforated grid block 

Simulator well model: 

► The well is defined by a list of 
connections 

► There is not notion of trajectory 

► Completions may be ordered by 
depth 



Calcul of simulator: 

► Gravity effect only 

► No friction or acceleration losses 

► Average density calculated for each 
segment 


BHP reference depth < 

► 




< 

> 



rfHh.. ... 

j 

k 




m 
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Wells representation: wells' Geometry 



Pressure at 



Average 

WBP9 

pressure 


DRAINAGE AREA 


Well cell 
pressure 
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Connection factor 


Definition 

► The production or injection rate of phase p for well connection k is: 


kr(S , .) 

n — CP ( D _ n \ 

Hp,k ~ K - Jjr well,k \^b,k r well,kJ 


► CF well is the Connection Factor, the well connection transmissivity or the 
numerical productivity index (numerical PI). 

► NOTE: 

• Pressures in the equation must be at the same reference depth 

• for a gas well, pseudo-pressure and non-Darcy flow may be considered, resulting in a 
different equation 
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Numerical PI 


Connection factor (numerical PI) 
Numerical pressure-rate relationship 



► For the phase p, the rate between two grid blocks 'i 1 and 'i+1 1 : 



► qp between the well and the block (b): 



► Pb and Sb are average values at the scale of the block, the pressure-rate 
relationship should take into account the physical gradients localised around 
the well and within the block volume: role of CF and kr(Sb) 
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Connection factor (numerical PI) 
Analytical pressure-rate relationship 



► Assuming radial flow, for steady-state and single phase: 


n — 

r 2 jr-K-h" 

(P(r 0 )-PJ 

Qp 

V? ) 

In (rjrj 



► If ro is such that P(ro) = Pb, the we derive the Connection Factor: 
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Connection factor (numerical PI) 
Analytical vs Numerical solution 




well half block 


CF ve ii — 


2k_ -Kh 

ln(r 0 /rJ + S 
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Natural PI vs numerical PI 



Numerical PI (CF): 

► Specific to the grid block containing 
the well 

► Units: m 3 .cP/(d.bar) 

► Independent of the phase produced 

► Constant in time 



B fx (Ln r b /r w + S + Cte ) 


Natural PI: 

► Related to the whole drainage area 

► Units: m 3 /(d.bar) 

► Depends on the phase produced 

► Varies with time 



B |x (Ln r e /r w + S + Cte ) 


Ln r,/r w + S + Cte 
nal Ln r b /r w +S+Cte 
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Connection factor (numerical PI) 
Peaceman's formula for ro 



► For a vertical well in the centre of a square grid block, Pb = P(ro) if: 

r 0 = 0.198 - Ax 

► Assumptions: 

• Vertical well 

• 2D 

• Pseudo-permanent, radial flow 

• Regular grid 

• Isotropic, homogeneous, non-dipping reservoir 

• Well far from other wells and reservoir boundaries 

• Well axis and axis aligned 


► If Dx = Dy, even considering anisotropy ratio Ky/Kx: 


r 0 =0.198 - Ax 
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Connection factor (numerical PI) 
Peaceman's formula for ro 

► Isotropic Reservoir - Square grid 
r w,b = 0,20 Ax 


► Isotropic Reservoir - Rectangular grid 

r »,b =0.I4 % /Ax 2 + Ay 2 

► Anisotropic Reservoir - Rectangular grid 


r - A-wVta k *)" Ax ' + ( k , k v)' : Ay ; 
w ’ b " a ‘ 8 (k y k,r + ( kx ky r 
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Numerical PI (CF) 
Importance 


► In History Match: calculation of ro in order to use Pbottom observations 

► From History Match to Forecast: 

• CF calibration by layers 

• Remember to update CF if grid cell permeability is changed (in H M) 

► In forecast: 

• use CF computed by the simulator from rw, S, block data 

• compute CF with SCFIEDULE 

• use available reservoir measures and adjust field PI 

• USE OF TOO HIGH VALUES MAY CAUSE NUMERICAL PROBLEMS!!! 
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Components of skin (remind) 


► The global skin factor represents various phenomena: 

• mechanical skin, Sm: local modification of permeability and effect of flow through the 
perforations. 

• geometric skin, Sg = Sd + Spp: effects of well geometry (deviation and partial 
perforation). 

• turbulence effect or quadratic skin, Sq: deviation from Darcy's law, representative of the 
turbulence effect. It mainly concerns gas wells but can be considered in the case of oil 
wells flowing at high velocity. 


s = s pp+ ^{s m+ s g+ s q ) 

n p 
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Bottom - surface relations 



P. 


CA 

P 

Y 

sep. 

j 

f Surface 



choke 



— # - Gas 

Gas metering 
P = separation pressure 
Oil metering 

— • oil 


Safety valve 
- Fluid produced 
Down hole choke 


PWF Well inflow 
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Inflow & outflow calculations 


► The total pressure drop has several components: 

• Pressure drop between the reservoir and bottom hole: "inflow" (1) 

• Well bore pressure head between connection and datum of VFP table (2) 

• Head losses in tubing: "outflow" (3) 
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Inflow & outflow calculations 



► Inflow performance relationship (IPR): section 1 

KV 

Q = CF—(P r -P wf ) 

M 

► Wellbore pressure head: section 2 

• Vertical well: hydrostatic head rgDz between the well's datum "connection i" and the 
VFP table datum depth 

• horizontal well : hydrostatic head rgDz + friction 

► Pressure drop in tubing: section 3 

• use of Vertical Flow Performance (VFP) tables provided in the data set 

• BHP = f (Q, THP, WCT, GOR, Qlift) 

► It is important to define accurately VFP tables when inflow drawdown is low 
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Outflow performance, VLP 



Pressure loss 
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Inflow & outflow 
Well deliverability 




Wells representation well control (gas lift) 




Liquid flow (nrP/d) 
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Wells activation - Bottom hole pumping 
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Inflow - outflow performance 



parameter 

IPR 

VLP 

Pres 

t 

+ 

no effect 

Pup 

T 

no effect 

- 

WCT 

T 

+ 

- 

GOR 

t 

- 

+ 

^tubing 

t 

no effect 

+ 

s 

t 

- 

no effect 

kro 

t 

+ 

no effect 



Positive effect on the VLP curve: 



there is an optimum value for GOR: a compromise between column 

weight and turbulence!!! 
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Outflow calculations 



► VFP tables are generally generated by well department 

► It is recommended 

• To use the PVT corresponding to produced fluids 

• To generate one table for each value of tubing diameter 

• To group wells according to their deviation, completion and produced fluid 

► The parameters must cover the range of values which will be reached during 
the forecast simulations: 

► NO EXTRAPOLATIONS ! 

► If extrapolation, some BHP may correspond to a combination of parameters 
physically impossible. 

• Values higher than the max reservoir pressure should be given to ensure that the well 
flowing conditions will not enter the forbidden region 




Matching of well performance 


► Well performance parameters are: 

• not used during history match (Wells are controlled in flow rate) 

• used during production forecasts (Well flow rates are directly linked to well 
performance) 

► During history match 

• It is pointless to try to match BHP and THP as long as reservoir pressure is not matched 

• It is highly recommended to match BHP and THP as soon as reservoir pressure is 
matched by introducing the VFP tables and adjusting PI. 

► At the beginning of production forecasts 

• BHP and THP must be matched per well to ensure the continuity from boundary 
conditions in flow (HM) to BC in pressure (forecast) 
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Well's constraints 

Production Primary and Secondary Constraints 



► Primary controls 

• Target rate of a principal phase 

• Fixed pressure: BHP or THP 

► Secondary controls 

• Maximum rates of one or more phases 

• Maximum ratios (GOR, WCT, WGR) 

• Limiting pressures 

• Limiting DP 

► Simulator will operate the well under primary control unless one of the 
secondary constraints is offended 

► Controls might be reset at any time during the simulation 
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Well's constraints 

Production Primary and Secondary Constraints 



► In practice there is no real distribution between primary and secondary 
constraints 

• Simulator will always honour the most constraining one 

► Example: if a well can produce 900 m3/d at BHP=30 bar 

• Target rate 1000 m3/d subject to minimum BHP of 30 bar: 

• Well produces 900 m3/d at BHP = 30 bar 

• Target BHP = 30 bar subject to maximum rate of 1000 m3/d 

• Well produces 900 m3/d at BHP = 30 bar 
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Well controls 



► Well controls 

• Control data for production and injection wells 

• Economic limit data for production and injection wells 

• Well efficiency factors (for downtime) 

• Maximum drawdown per producer 

► VLP tables for production and injection wells 
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Well controls in history match 



► In history Matching, observed average rates are known; controls are simple. 

► VFP tables are introduced at the end of history matching process to ensure 
the continuity between matching runs (set measured Q) and prediction runs 
(limit THP). 
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Simulator computations 



► Well dispatch of group constraints: 

• evenly 

• as specified by the user 

• according to well potentials (PI) 

► Rate calculation for all phases using 

► Pressure imposed: 

► From Phead to Pbottom with pressure loss tables 

► Search of stable functioning point: iterative process 
¥ Rate calculation per phase and per perforated cell 
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Group of well constraints 


► Hierarchy: Field, Platforms, Groups, Wells 

► Constraints: 

• On rates for main phase (economical, contracts, capacity) 

• On rates for secondary phase (capacity of flaring, water treatment) and gas lift rate 

• X-dependency: Qproduction = f(Qinjection), Pressure maintenance 

► Monitoring: 

• Opening of waiting wells 

• Shut-in of most water flooded well 
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Field and numerical PI 
Exercise 

► Calculate the numerical PI and the field PI for a monophasic vertical oil well 
from the following data: 

• Qo = 1000 Sm3/d 

• KH = 1000 mD.m and S = 0.5 from well test interpretation 

• rd = 1000 m 

• rw = 0.1m 

• po =0.7 cP and Bo = 1.2 

• The well is in a square grid of side 100 m and the formation thickness is 10 m 

• Interference test show a permeability anisotropy kx/ky = 0.5 

► If the bottom hole flowing pressure is 100 bar, calculate: 

• Reservoir pressure at the drainage radius 

• The well grid block pressure (Pb) in the reservoir model 
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Keypoints to keep in mind 

CXu 

\ 

► Main constraints for the simulator: 

• Pressure target: 

• Maximum or imposed 

• WHP/BHP 

• Minimum for producer - maximum for injector 

• Maximum for some layers 

• Rate target: 

• Surface rates: Qo, Qg, Ql 

• Bottom: Qt (history match) 

► Secondary constraints for the simulator: 

• Water cut 

• Gas oil ratio 

► Monitoring: actions when a constraint is offended: workovers? 
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Courlis oil field was discovered in 1984 under 80-90 meter water at about 30 km from the coast line. 

The first exploratory well, Courlis Marine 1 (CLM1), targeted a deep structure. During drilling, direct 
fluorescent and samples showed bitumen and residual oil before the targeted structure was reached. They 
were located in Cenomanian sands. Small gas indices were also found in Turonian carbonates and Senonian 
clay-rich sandstones. 

A few kilometers from CLM1, a shallow anticlinal structure showed bright spot in the Senonian, a possible 
flat spot in the Cenomanian and a fuzzy zone near the top of the structure. Thus, it was decided to begin a 
high resolution seismic campaign and drill a second well, CLM2. This well showed oil reservoir in Cenomanian 
and Turonian series and gas in the Senomian which was very shallow. Furthermore, the water-oil contact was 
found in the Cenomanian at depth corresponding to the flat spot. 

The structure area (18 km2), the good well productivity from CLM2 and the good reservoir characteristics for 
Cenomanian emphasized the interest of the discovery. 

During the first quarter of 1984, 4 evaluation wells (CLM3, CLM4, CLM5, CLM6) were drilled to better define 
the structure boundaries and check that flat spot matches the water oil contact. From these wells, it follows 
that the flat spot was due to lithology changes. However, the structure size was confirmed. Gas was detected 
in the Senonian, oil in Turonian and Cenomanian reservoirs. 

This exercise will only deal with the Cenomanian reservoir. 
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Cenomanian reservoir description 


At the top of the Cenomanian, the structure is an anticlinal along a NW-SE axis with a narrow top graben. The 
reservoir top is at 550 mSL. CLM6 well was almost at the structure top. A slightly higher zone may exist North 
West of CLM6. 

The Cenomanian reservoir has three level (Cl to C3). We are mainly interested in the top of Cenomanian 
series (Cl level) which is made of clean unconsolidated sands with fine to coarse grain inter-bedded with silt 
and shale. The top of Cenomanian reservoir (Cl level) is located at a depth of 545 mSS and the OWC is at 645 
mSS. The thickness of Cl level varies between 10 and 45 m with a maximum near the centre of the structure. 

At varying depth from the top of the reservoir, the C2 level is made of thigh dolomitic sandstone with fine to 
medium grain sandstone both inter-bedded with shale. The thickness of this level varies but it remains quite 
compact. The C2 reservoir is quite thick at the centre (over 90 m). Its thickness decreases down flank to zero. 
It created the seismic amplitude variations at the Cenomanian top. This layer presents good quality 
sandstone but they are not continuous (see CLM2 log). The transition between this layer and sands is not 
well defined in CLM6 well. The C2 level is essentially oil bearing and core plugs (CLM3, CLM4, CLM5, CLM6 
wells) show that this is an a priori tight level. 

Under this C2 level, there are thick aquifer sands (C3 level). The aquifer reservoir C3 is about 250 meter 
thick. However the Cl reservoir is fully isolated from C3 aquifer due to the a priori tight dolomitic sandstone 
in C2 reservoir. 
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Isobathes top cenomanien 






Reservoir characteristics: 

Cenomanian petrophysical characteristics 



► Using log interpretations for CLM2, CLM3, CLM4, CLM5, CLM6, the geological 
study could establish correlation between these wells and a geological model 
was build for the field. Log interpretations yielded porosity and initial water 
saturation as showed in the following table: 


WELL 

Top Cl 
reservoir 
fmSS) 

Top C2 
reservoir 
(mSS) 

h, 

(m) 

K 

(m) 

<i» 

(%) 

Swi 

(%) 

CLM2 

588 

615.5 

27.5 

20 

29 

11 

CLM3 

629 

648 

19 

15.5 

28 

11 

CLM4 

627 

637 

10 

8 

27 

11 

CLM5 

629 

637 

8 

6 

26 

11 

CLM6 

555 

596.5 

41.5 

30 

28 

11 
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Fluid properties 



► Oil 


► Water 


Property 


Reservoir temperature 

46° C 

Initial reservoir pressure 

59 bar 

Saturation pressure 

57 bar 

Oil formation volume factor 

1.105 v/v 

Solution gas @ Psat 

32 v/v 

Viscosity @ Psat 

10.6 cP 

Compressibility @ Pb 

0.000109 bar 1 

Stock tank oil specific gravity 

0.896 

Property 



Water formation volume factor 
Density @ 25° C 
Viscosity 
Compressibility 


1 .01 v/v 
1.07 
0.7 cP 

0.000035 bar 1 
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Initial conditions 



► The initial reservoir temperature is 48°C and its initial pressure is 59.1 bar at 
the reference depth of 583 mSS. Therefore, the reservoir is slightly under 
saturated and there is no clear evidence of initial gas cap. Furthermore, no 
variation of bubble point versus depth was noticed. The initial water oil 
contact is at 645 mSS. 

► Water-Oil and Gas-Oil relative permeability curves were estimated from 
similar reservoirs. 

► Since the rock permeability is quite high, a hysterisis phenomenon on the 
capillary pressures was considered. The rock compressibility was estimated to 
be 5 10-5 bar-1. 
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W/0 relative permeabilities 


Irreducible water saturation 

Swc 

= 11% 

Residual oil saturation to water 

Sorw 

= 23% 

Relative oil permeability at Swckro 

= 1 


Relative water permeability at 1-Sorw 

krw 

= 0.32 
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G/0 Relative Permeabilities 


Critical gas saturation 

Sgc 

= 2% 

Residual oil saturation to gas 

Sorg 

= 23% 

Relative oil permeability at Sgc kro 

= 1 


Relative gas permeability at 1-Sorg-Swc 

krg 

= 0.3 
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Numerical model 


► In order to propose the development scenario, a 3D, 3-phase, black oil 
numerical model was built for the simulator. The horizontal grid which is 
oriented North East - South West is 31 by 12 in X and Y direction respectively. 
Each cell is 250 m by 250 m. Since no zonation was established in Cl reservoir, 
this level was divided in 3 equal thickness layers. 

► Petrophysical characteristics and fluid data for the model are obtained from 
the previous information. They are already provided in the data file. 

► Then you should propose development scheme through at least a primary 
production to satisfy an annual production requirement of 2% Initial Oil In 
Place as long as possible. However, during early development stages, 
production should be carefully planned to maintain reservoir balance. A 15 
year production forecast is desired. 
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Grid-top reservoir 




!FP Training 


467 


Contact 


Simulation grid (oil saturation) 




OMSfll 



OOOOM 


0 21250 


044500 


0 60750 


0,89000 


3 layers, nx = 31, ny = 12 => 1100 active cells 
No aquifer activity 
Net to gross ratio: 0.75 
Average porosity: 27% 

Average Swi : 11% 

Horizontal permeability: 1400mD 

WOC @ -645 m 

Top reservoir: -561 m 
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Project objectives 



► Initialize the numerical model. 

► Compute the initial fluids in place in Cl reservoir. 

► 

► Several development schemes must be investigated for Courlis Marine Oil 
Field: 

• Produce the field through natural depletion finding the optimum number of wells for 
this case. Make a sensitive study: Change Sgc to 10% and draw some conclusions. 

• Investigate secondary oil recovery using water injection: which is the main problem 
encountered. 


IFP Training 


469 


Project constraints 


► Wells are drilled from a maximum of 3 platforms A, B, C. The average well diameter is 81/2". Each platform 
has a maximum of 10 wells. Wells have a 31/2" tubing. 

► For production wells, a 10 bar initial drawn down at sand face is the minimum needed to operate the well. 

No vertical flow performance tables are provided. 

► Horizontal and vertical wells may be considered for production. 

► Well drilling and completion take about 20 days. Two rigs are available during the initial field development. 
Wells are put on-line by group (4 wells per group). 

► The bottom hole pressure is 25 bar for the producer during the natural depletion. 

► The maximum admissible production GOR is 500 m3/m3 and the maximum allowable water cut is 95%. The 
minimum economical rate to operate a well is 5 m3/d. 

► It may be necessary to optimize the perforation intervals to enhance recovery. 

► Each well may produce up to 300 m3/d oil for a vertical well and up to 1000 m3/d for a horizontal one. 

► For secondary recovery, sea water may be used since no compatibility problem exists. The fracturing 
pressure is 100 bar for this reservoir. The injection flow rate should be less than 400 m3/d and selected to 
optimize recovery. 

► A minimum spacing of 500 m should be used between injection and production wells to avoid early water 
break through. 
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Reservoir simulation workshop 
History match 
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EXPLORATION APPRAISAL/ 
DELINEATION 


FIELD CONSTRUCTION 


History match 


Objective 

► Validate (Calibrate) the reservoir characterization by comparing performance 
of the model with historical performance (rates, pressures, saturations) 

Principle 

► Reproduce with the model the measured evolutions of pressure, BSW and 
GOR by well, by zone or for the entire field 

• If discrepancies appear during this comparison, 

• Modify the characterization in a rational way 

• Change uncertain parameters first (e.g. aquifer) 

• If necessary, change the characterization 

• Thus at the end of history match, the characterization should be improved 

► The history matched model can be used for prediction with more confidence 

• Difficulties: uncertainties on fault and flow barriers network 

• Rule of thumb: predictions are reliable on a period twice the production period 
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History matching: main issues 


► Observed flows are imposed on wells during the history period 

► One expect to reproduce: 

• Pressure evolution 

• WCT and GOR 

• Gas or water breakthrough 

• WOC and GOC contacts 

► This is not a simple work as: 

• Many data are unknown (little or no information is available at distance from wells) 

• It is not obvious to detect the most influent parameter (all parameter act together) 

• Some artefacts must be corrected (grid size, grid orientation...) 

► It is possible to distinguish between two main types of problems: 

• Pressure match 

• Saturation match 
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History match 
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History matching: methodology overview 



► 4 main steps 

• Step 1: Identification of available data that have to be matched 

• Adapt data to grid size 

• Step 2: Data Analysis 

• Identification of main uncertainty in the Geomodel 

• Step 3: Selection of matching parameters 

• Identification of probable range for each matching parameter 

• Step 4: Modification of matching parameters 

• Trial & error process 

► Reminder 

• G&G must work hard to help the reservoir engineer to maintain the consistency of the 
geological model 

• It is better to have rough, consistent matching than matching which is accurate but 
destroys the model 
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History matching: schematic workflow 




NEW GEOMODEL 


YES 
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History match: data to match 
Production Data - Well rates 



► Production data 

• Field production (Oil, Gas and Water) 

• Allocated production per well (Oil, Gas and Water) 

• Allocated production per set of perforations (if PLT available) 

► Pressure data 

• Type of pressure measurement 

• Timing and frequency of measurement 
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History match : data to match 
Well rates 

► Determine the accuracy and frequency of measurements 

• Which kind of separator (Gas-Liquid, 3 phase separator) 

• How is measured WOR, GOR if no 3 phase separator? 

• Frequency of measurements (weekly, monthly?) 

• Activity factor (% of activity of the wells) 

► Allocation of rates to perforated interval 

• Completion and well status (casing, cement) to be known 

• Production logging is the best tool to allocate rate 

• If no production logging, estimate the accuracy of allocation to intervals 

• When possible, draw the maps of injected fluid breakthrough for each interval 
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History match: data to match 
Pressures - Saturations 

► Pressures 

• Shut-in pressures of wells 

• Look at pressure curves to estimate the pressure in the cell 

• RFT in wells drilled after the start of production 

• Differential depletion by interval 

• Communication through faults 

• Observation wells 

• Shut-in wells should be changed to observation wells 

• WHP are usually not used 

• flow in tubing difficult to match exactly, interference with surface flow lines 

► Saturations 

• Cased Hole logging, PLT, coring 
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History matching: some advice 



► Start with a geological model 

► Determine the level of uncertainty related to the different parameters 

► At each step during the match process, note the influential parameters in 
order to limit the number of runs 

► Use analytical model to evaluate the effect of each parameter and determine 
the most influent 

► Take care with gridding effects 

► Run sensitive cases for the influent parameters 

► Adjust the primitive analysis when the match remains poor 
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Strategy 




Strategy 

Field and well basis match 

► Field basis match: 

• Faults 

• Aquifer 

• Global permeability scaling 

• Vertical Transmissivity 

► Well by well match: 

• Local Transmissivity X, Y, Z 

• Relative permeability endpoint scaling (Swi, Sor) 

• Local PI and skin 
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History matching strategy 


First phase: pressure match 

► Look at total fluid production and average pressure (kind of material balance) 

• Of the whole field 

• Of parts of the field (compartments, or zones determined from the geological analysis, 
layers) 

► Change first the more uncertain parameters by zone 

• Aquifer transmissibility (kh), storativity (kh ct) 

• Reservoir permeability 

• Multiplying factor to reproduce pressure gradients 

• Vertical connections to account for pressure discrepancies between layers (RFT useful) 

• Connections through faults to account for different pressure regimes 
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History matching strategy 



Second phase: Fluid distribution (saturation match) 

► By zones (compartments, or zones determined from the geological analysis, 
layers) look at contacts movements, WOR and GOR of wells. Try to match 
fluid BT (break-through), fluid produced volumes 

► Adjust first the permeability distribution 

• Vertical distribution by layer and connection between layers (vertical permeability) 

• Areal distribution of permeability (barriers, high permeability zones, sealing or 
conductive faults) 

► Change Kr only if changes in permeability distribution cannot achieve a 
satisfactory match 

• First check if initial water saturation is correctly represented 

• Look at the possibility of hysteresis 
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First step - General field match - Run 1 



Wrong initial 
Pressure drop 



7 rung 491 


Next step - General field match - Run 2 
Aquifer potential increase 



Good Initial 
Pressure drop 



0000 3 


0— 
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Final step - General field match - Run 3 



Good match of 
pressure profile 


Oil production 
match 


Oil permeability increase 



100^ 

80 3 
60 ^= 
40 ^ 

20 3 
0 


Water Celt (%) 

Moderate Water Prod. 
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History matching strategy 



► Final step: well by well match 

• Once the global and zonal matches are correct, look at each well 

• Check if the cell size is not the cause of an incorrect match 

• Check if coning can be suspected and is not taken into account by specific well Kr 
functions 

• Check if the discrepancy does not reveal a completion problem (cement or casing leak, 
fluid entry from another interval) 

• Corrections should remain in the vicinity of the well 
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History matching strategy 
Final step: well by well match 




History matching: some advices 



► Flow directions are not correct if pressure is not matched 

• Don't try to match saturations if you are not matched in pressure 

► Early well behaviour correspond to area close to the wells 

• Concentrate on well data to match early production times 

► Late well behaviour correspond to area far from the wells 

• Don't limit your analysis close to the wells to match late production times 

► Modification of matching parameters 

• Try to anticipate model reactions by using simple calculations 

• Don't introduce new parameters without a look back to G&G 
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Pressure match: material balance 


► Objective 

• Get a correct evolution with time of the average reservoir pressure 

► Main parameters 

• Volumes Originally in Place 

• Aquifer size & water influx 

• Pore & Fluid Compressibility 

► Important notice 

• The material balance should address the whole reservoir voidage (no material balance 
per fluid in surface conditions) 

• It is useful to get an energy balance to have an estimation of the importance of each 
individual production mechanisms (pore volume contraction, fluid expansion, water 
sweep, gas sweep...) 
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Pressure match: material balance 



► Pressure match is an adjustment of the reservoir energy balance between: 

• Volumes originally in place 

• Aquifer activity 

• Pore and fluid compressibility 

► Reservoir voidage calculation 

Qres = Qo-Bo + Qg Bg -Rs-Qo-Bg + Qw- Bw 


► Important notice 

• Reservoir voidage has to be calculated and is depending on the reservoir pressure 

• ECLIPSE keyword for reservoir voidage is RESV 
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Pressure match: Darcy's law 


► Objective 

• Get a correct geometry of the flow lines and pressure drop along flow lines 

► Main parameters 

• One phase flow 

• Transmissivity distribution 

• Multi phase flow 

• Transmissivity distribution 

• Transfer functions (relative permeability & capillary pressure) 




Pressure match 
Diffusivity equation 



Objective: Get a correct evolution of reservoir pressure versus time and space. 


► Diffusivity equation: 

d 2 P d 2 P d 2 (P- p-Q'Z) _ K dP 
dx 2 dy 2 dz 2 </>•) u-c dt 


► Main parameters: 

• hydraulic diffusivity, K/(4>.|j..c) 

• permeability, K 

• fluid viscosity, \x 

• porosity, § 

• total compressibility, c 
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Pressure match: wells' representation (1) 



► Objective 

• Get a correct relationship between flow rate, reservoir pressure and bottom hole 
flowing pressure 

► Main parameters 

• Numerical productivity index 

• Drainage area properties 

• Transmissivity distribution 

• Transfer functions (relative permeability & capillary pressure) 
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Pressure match: wells' representation (2) 



Pressure at 

drainage 

radius? 


Bottom hole 

flowing 

pressure 



Average 

WBP9 

pressure 


Well cell 
pressure 


AREA REPRESENTED BY THE NUMERICAL PI 
DRAINAGE AREA 
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History match strategy 



► First phase: pressure match 

• look at total fluid production and average pressure (kind of material balance) 

• of the whole field 

• of parts of the field (compartments, or zones determined from the geological analysis, layers) 

• change first the more uncertain parameters by zone 

• aquifer transmissibility (kh) , storativity (kh ct) 

• reservoir permeability 

• multiplying factor to reproduce pressure gradients 

• vertical connections to account for pressure discrepancies between layers (RFT useful) 

• connections through faults to account for different pressure regimes 
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History match strategy 



► Pressure match parameters are: 

• Depending on hydraulic diffusivity K = k / <Dpc 

• Mobility k/p 

• Porosity <1> 

• Compressibility c 

• External gridding limit conditions (aquifers, flow transfer ...) 

• Pstatic evolution with withdrawal -> adjustment of O and c (AP = 1/c * AV/V) 

• Increase of O or c => AP decrease) 

• Pstatic - Pflowing evolution -> adjustment of k/p (Darcy) 

• Increase of k/p => AP decrease 
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Pressure match 
Aquifer activity 



► A preliminary study with a MB software is necessary to run the aquifer 
match. The expected result is the aquifer volume plus its permeability. 


► The aquifer volume is to be reproduced in the reservoir model with aquifer 
cells or analytical functions. 


► Aquifer activity needs to be adjusted in order to reproduce the field observed 
reservoir pressure history. 


► The reservoir model production history is run with all the producing wells 
governed by the "reservoir voidage" option. 
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Pressure match 
Fault modelling 




► Fault modelling: Only faults with influence in the zone of interest are 
modeled. 

• definition 

• transmissivity 
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Pressure match: wells' representation 
Numerical PI - Analytical vs numerical solution 




well half block 



2 n-K-h 
ln(r 0 /r w ) + S 
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Flow rates match 



► During the pressure match procedure, production rates are not honoured. 

► Simulation is now run by setting the oil rate for producers so bottom hole 
flowing pressure, gas and water rates are calculated by the simulator. 

► The phase rate matching consists of adjusting the calculated GOR and WCT to 
the field measured values. 

► To honour the relationships between reservoir pressure, BHFP and phase 
rates, PI need to be adjusted. This is accomplished by applying multiplication 
factors to the well perforation connection values: MULTPI. 


/PPTraining | 509 

Flow rates match 

► Main parameters: 

• Permeability -Transmissivity. 

• Permeability barriers (i.e. faults) 

• Relative permeability: shape and endpoints. 
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Flow rates match 
Relative permeabilities 



k\ 

- end-point relative 


permeability for oil 

k\ 

- end-point relative 


permeability for water 

Kc 

- connate water saturation 

9 

o res 

- residual oil saturation 



Anatoly Zolotukhin, 
Statoil, 2004 


!PP Training 


511 


Flow rates match 
Remind: Endscale option 



► Objectives 

► Modify relative permeabilty tables in an easy way, kr tables are normalized 
and remain always the same, only the end-points are changed and kr curves 
are then recalculated. 


► It's a useful option in History Match simulations. 
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Flow rates match 



Remind: Endscale option 

KEYWORDS 

► KRO: Maximum oil relative permeability 

► KRORW: Oil relative permeability at 
critical water saturation Swcr 

► KRWR: Water relative permeability at 
residual oil saturation (1-Sowcr) 

► KRW: Maximum water relative 
permeability 

► SWL: Connate water saturation 

► SWCR: Critical water saturation 

► SOWCR: Residual oil saturation 

► SWU: Maximum water saturation 



Flow rates match 



Relative phase permeabilities 

► Water cut (WCT): water - oil ratio 

► Water breakthrough (WBT): water production starts - associated to Swi 


w v FWCTH vs TIME (Slnn_run) 
”FWCT vs. TIME <Slm_run) 
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Flow rates match 
Water breakthrough 

► Matching breakthrough times is a difficult task. 



► Breakthrough times are sensitive to truncation errors (numerical dispersion) 
and the accurate matching requires finer grid than normally necessary. 


► Using a LGR is a possibility to the use of pseudo-relative permeabilities can 
help. 

► An unsuccessful attempt for a match indicates that some of the basic 
assumptions of the model (geology, structure, volumes, extensions, PVT 
behaviour, energy balance between initial hydrocarbon in place and aquifer 
activity) may have to be revised. 
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Flow rates match 
Water breakthrough 



► The gridding techniques include local gridding (LGRs) for the creation of small 
cells around wells for improved resolution, useful to match the water 
breakthrough and water cut when conning effects are present. 
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Flow rates match 
Water breakthrough 




Cumulated Water (bubble diagram from OFM) 
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Uncertainty contributions 



► Seismic interpretation 

► Geological interpretation 

► Well production allocation 

► Observed data (bottom hole pressures, well logs, etc...) 

► Fluid models 

► Simulation grid accuracy (e.g. fault juxtapositions), orientation, etc. 

► General poor control outside the main reservoir structure (aquifer support 
etc...) 
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Key points to keep in mind 


CXu 


B 

► Keep a coherent geological model 

► Do not match one parameter while ignoring the others 

► Do not discard data which do not match!!! 

► Avoid multiplying the Kr curves 

► Do not look for a perfect match 

Be very patient! 

V 
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Reservoir simulation workshop 
Production forecast 



Production forecast 


► Introduction 522 

► Well controls 524 

► Inflow & outflow calculations 528 

► Matching of well performance 538 
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Introduction 



► Good Production forecasts will depend on 

• Good history match 

• Good integration of well performance and production constraints 

► Well performance 

• Inflow curves 

• Productivity index or complete curve 

• Outflow curves 

• VFP tables 

► Production constraints 

• Constraints related to flow 

• Maximum water-cut or maximum GOR per well 

• Maximum water production or gas production for a group of wells 

• Constraints related to pressure 

• Minimum bottom hole flowing pressure 

• Minimum well head pressure 
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Well Controls: principal uses in forecast mode 



► Main controls 

• Imposed Flow rate 

• Oil, Water, Gas, Liquid, Reservoir voidage 

• Imposed pressure 

• Tubing Head Flowing pressure 

• Bottom hole flowing pressure 

► Secondary controls 

• Flow rates 

• Economic limits 

• Rates upper limit per phase 

• Maximum ratios (Wcut, GOR, WGR) 

• Pressures 

• BHP and THP limit (lower for a producer, upper for an injector) 

• Maximum drawdown 

► VFP tables are used to relate bottom hole to well head pressures 
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Well controls: ECLIPSE keywords 


► Well controls 

• WCONPROD: Control data for a producer 

• WECON: Economic limit data for a producer 

• WELLDRAW: Set maximum drawdown per producer 

• WCONINJE: Control data for an injector 

• WECONINJ: Economic limit for an injector 

► VLP table 

• VFPPROD: VFP tables for a producer 

• VFPINJ: VFP tables for an injector 

► Dimensions in RUNSPEC 

• VFPPDIMS: Production wells VFP table dimensions 

• VFPIDIMS: Injection wells VFP table dimensions 
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Inflow & outflow 
calculations 



Inflow & Outflow calculations 



► Inflow 

• Pressure drop between cell and bottom hole (section 1) 

► Outflow 

• Pressure head between connection and datum of VFP table (section 2) 

• Head losses in tubing (section 3) 



Deviated well 


Horizontal well 


Inflow & Outflow calculations 


► Inflow performance relationship (section 1) 

kr 

Q = CF x — x (Pr — Pf ) 

[A. 

¥ Wellbore pressure head (section 2) 

• Vertical well: hydrostatic head between the connection i and the well's datum + 
(correction up to VFP datum depth) 

• Horizontal well: hydrostatic head + friction if WFRIC or MSW 

► Pressure drop in tubing (section 3) 

• Use of Vertical Flow Performance (VFP) tables provided in the data set 

• BHP = f (Q, THP, Wcut, GOR, Qlift) 

► Reminder 

• It is important to define accurately the VFP tables when the inflow drawdown is low 
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Outflow calculations: VFP tables generation 



► VFP tables are generally generated by well department 

► It is recommended 

• To use the PVT corresponding to produced fluids 

• To generate one table for each value of tubing diameter 

• To group wells according to their deviation, completion and produced fluid 

► The parameters must cover the range of values which will be reached during 
the forecast simulations: NO EXTRAPOLATIONS! 

¥ Some BHP may correspond to a combination of parameters physically 
impossible. 

• Values higher than the max reservoir pressure should be given to ensure that the well 
flowing conditions will not enter the forbidden region 
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Matching of well 
performance 


Matching of well performance 



► Well performance parameters are: 

• Not used during history match (Wells are controlled in flow rate) 

• Used during production forecasts (Well flow rates are directly linked to well 
performance) 

► During history match 

• It is pointless to try to match BHP and THP as long as reservoir pressure is not matched 

• It is highly recommended to match BHP and THP as soon as reservoir pressure is 
matched by introducing the VFP tables and adjusting PI 

► At the beginning of production forecasts 

• BHP and THP must be matched per well to ensure the continuity from boundary 
conditions in flow (HM) to boundary conditions in pressure (forecast) 
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Matching of well performance: BHP per well 


► BHP per well can be modified by 

• Tuning of Connection factors (CF) 

• Introduction of production uptimes 

• Use observed daily Q. + WEFAC instead of monthly averaged Q 

• Use of WPIMULT or modification of the skin 

► Reminder: Tuning CF per layer may be necessary for PLT match 

• The information provided by the PLT (allocation of flow per layer) must be taken into 
account since the match in P (K per layer) 

• It is recommended to adjust the skin per (group of) layers 

• Use instantaneous rates (not monthly averaged) and be sure that everything (P, 
measure conditions) is consistent 
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Matching of well performance: THP per well 


► BHP per well can be modified by 

• Modification of the pressure head between connection and datum of VFP table 

• Modification of the VFP table 

► Eclipse solves for the intersection of: 

• The linear portion of the VFP table at the flowing conditions (GOR and Wcut) 

• The local IPR: BHP = Pr - (Q/IP) 

► This intersection may be in the unstable part of the VFP curve; in this case the 
results are very sensitive to the slope of IPR 

► It is recommended to discretise each parameter in VFP curves with enough 
values to represent correctly the working point 
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Thank you for your attention! 


